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1 Introduction 

The minimal renormalizable supersymmetric SU(5) model [1], with just 3 pairs of 10 © 5 
fermion representations and an adjoint 24 plus a 5 © 5 pair in the Higgs sector, is the 
simplest supersymmetric Grand Unihed extension of the Standard Model. It is therefore 
particularly important to test this model in detail, and possibly to rule it out. Although the 
choice of the gauge group, supersymmetry and minimality do not need a special motivation, 
it is more difficult to justify the absence of non-renormalizable terms in the superpotential. 
Experimental evidence tells us that some of these terms must be strongly suppressed. For 
instance, the superpotential operators QiQiQ 2 Lk/Mpianck (k = 1, 2, 3) induce proton decay 
at an unacceptable rate unless they come with coefficients smaller than about 10“^. Their 
smallness will lead us to assume in this paper that for some (to us unknown) reason all 
non-renormalizable operators can be neglected, and to adopt the minimal renormalizable 
supersymmetric SU(5) model as a benchmark^. 

It has been shown long ago [7] that the region of the parameter space of the minimal 
renormalizable supersymmetric SU(5) model corresponding to TeV-scale soft terms is ex¬ 
cluded^. The reason put forward was the incompatibility between the colour triplet mass 
constraints associated with gauge coupling unihcation on the one hand, and with proton 
decay on the other hand. This conclusion relies however on the assumption of a relatively 
light superpartner spectrum (although masses as large as 10 TeV for the hrst two genera¬ 
tions of sfermions were considered in Ref. [7]). The purpose of this paper is to ascertain 
whether it can be extended to the region of larger superpartner masses. 

In order to answer this question, several constraints have to be imposed on the model: 
precise gauge coupling unihcation, correct predictions for charged fermion masses, the Higgs 
mass constraint, the experimental bound on the proton lifetime, and hnally the experimen¬ 
tal bounds from Favour physics and from direct searches for supersymmetric particles. We 
also require perturbativity of the model. The observed down-type quark and charged lepton 

Another option is to give up renormalizability and to assume that the non-renormalizable operators 
giving rise to fast proton decay are suppressed, while harmless higher-dimensional couplings can be sizable. 
Under this assumption, it is possible to avoid fast proton decay from heavy colour triplet exchange [2, 3, 4, 5] 
and to correct the (phenomenologically inaccurate) SU(5) relations between charged lepton and down-type 
quark masses [6] (for a review on these issues, see Ref. [4]). 

^In the case of decoupling scenario of heavy first two generations of sfermions considered in [7] proton 
decay can still be consistent with the experimental bounds providing the flavor sfermion sector gets a very 
specific form [3]. 
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masses are accounted for by generation-dependent supersymmetric threshold corrections 
with large A-terms [8, 9, 10, 11, 12, 13, 14], which in turn are bounded by considerations 
related to the stability of the electroweak vacuum [15, 16, 17]. The soft supersymmetry 
breaking terms are required to respect SU(5) invariance but are not assumed to be flavour 
universal as in the so-called mSUGRA model (as a matter of fact, generation-dependent 
A-terms are needed to correct the SU(5) fermion mass relations). To avoid potentially 
dangerous flavour effects, we will therefore take the GUT-scale sfermion soft terms to be 
aligned with fermion masses^ [18]. Neutrino masses can be generated either through bilin¬ 
ear R-parity violation"^ [19, 22, 23, 24] or by adding right-handed neutrinos to the model in 
order to implement the seesaw mechanism [25]. Finally, if the neutralino sector does not 
contain a suitable candidate for dark matter (either because the lightest neutralino is not 
the lightest supersymmetric particle, or because it is too heavy), this role may be played 
by the gravitino. 

Not surprisingly, the main constraint on the superpartner mass scale comes from pro¬ 
ton decay, which pushes the spectrum above the 10 TeV scale (the conflict between the 
proton decay and unihcation constraints pointed out in Ref. [7] is resolved by relaxing the 
upper bound on the superpartner masses). Perturbativity imposes an upper bound on the 
coloured triplet mass, which translates into an upper limit on superpartner masses once 
gauge coupling unification is imposed. The observed Higgs mass, together with vacuum 
metastability constraints associated with the stop A-term, also excludes large portions of 
the parameter space. A priori, there is no guarantee that the minimal renormalizable su¬ 
persymmetric SU(5) model will survive all these constraints, even if very large values of 
the soft terms are allowed. 

To give an idea of how the parameter space is restricted by the various phenomenological 
requirements, we show in Fig. 1 the approximate constraints in the (tan/?, rrisusy) plane 
obtained by making several simplifying assumptions. Namely, all sfermion masses are taken 
to be equal to nisusy at the low scale, as well as the /i parameter and rrin^ > 0 is 

assumed, and and the B parameter are computed from the electroweak symmetry 
breaking conditions), while the gaugino masses assume a common value M 1/2 = msusy at 
the GUT scale and are split by renormalization group running. Obviously, these inputs 
are not consistent with SU(5) symmetry of the soft terms at the GUT scale, but they 
make it possible to show several constraints in a single plot. Imposing SU(5) boundary 
conditions will signihcantly affect quantities such as the heavy colour triplet mass and more 
crucially the proton lifetime, making the investigation of the parameter space of the minimal 
renormalizable supersymmetric SU(5) model more involved than suggested by Fig. 1. As 
we are going to see, phenomenologically viable points typically feature superpartners in the 
0(10^ — 10^) TeV range, with values of tan/? between 2 and 5, but lighter spectra with 
supersymmetric particles as light as a few 10 TeV can also be found. 

In the process we have generalized the procedure of Refs. [26, 27] for deriving approx¬ 
imate semi-analytic solutions to the one-loop renormalization group equations (RGBs) for 
the MSSM soft terms. In this way we are able to write the low-energy soft terms as linear 

^Flavour violating soft terms are going to be generated from renormalization group running (due to the 
CKM matrix, and possibly to right-handed neutrino couplings and R-parity violating couplings), but this 
will lead only to small effects in flavour-changing processes, especially in view of the large superpartner 
masses. We shall therefore neglect these small RG effects. 

^The SU(5)-invariant bilinear R-parity violating operators also contain baryon number violating 

terms, which however are harmless by virtue of the doublet-triplet splitting if <C Mgut [19j 20]. Contrary 
to Ref. [21], we assume negligible trilinear R-parity violating couplings. 
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Figure 1: Approximate constraints on the minimal renormalizable supersymmetric SU(5) 
model in the (tan/3, plane. The region to the left of the red curve is excluded by the 

measured Higgs boson mass, while the area to the right of the purple curve is excluded by 
vacuum stability considerations (namely, the value of the stop A-term needed to reproduce 
the observed Higgs boson mass would render the electroweak vacuum unstable with a 
lifetime shorter than the age of the universe). The experimental lower bound on the proton 
lifetime is satished above the green line, and the perturbativity constraint tut < Mpianck/^^, 
where niT is the colour triplet mass, is satished below the black line. The green area is 
consistent with all four constraints. In this plot, the sfermion masses, the fi parameter 
and mHy, are all equal to msusy, as well as the common gaugino mass at the GUT scale. 
Since these values are not consistent with SU(5) boundary conditions, the hgure should be 
considered as illustrative only. 

or quadratic functions of the initial (GUT-scale) parameters, making it possible to explore 
the parameter space without having to solve the RGBs for each point. In practice, one just 
needs to solve numerically the RGBs for gauge and Yukawa couplings for each choice of 
tan/3 and rUsusy, the matching scale between the SM and the MSSM. The low-energy soft 
terms and their dependence on the other model parameters are then simply given by linear 
and quadratic algebraic equations. 

The paper is organized as follows. Section 2 introduces the minimal renormalizable 
supersymmetric SU(5) model. In Section 3, the running of the model parameters is dis¬ 
cussed, and the semi-analytical procedure used to solve the renormalization group equations 
(RGBs) is presented in Section 4 (more details can be found in Appendix A). Proton decay 
and the constraints associated with the metastability of the electroweak vacuum are ad¬ 
dressed in Sections 5 and 6, respectively, with technical details relegated to Appendices B, 
G and D. Finally, we present our results in Section 7. 
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2 The minimal renormalizable supersymmetric SU(5) 
model 

In this section, we briefly describe the minimal renormalizable supersymmetric SU(5) 
model [1] and present our notations. The Higgs sector includes the adjoint 24//, which 
spontaneously breaks the SU(5) gauge group to SU(3)xSU(2)xU(l), and a fundamental 
and anti-fundamental representations 5// and 5// containing the two light Higgs doublets 
responsible for electroweak symmetry breaking. The bn © Higgs fields also includes a 
heavy pair of colour triplet and antitriplet that mediate proton decay through d = 5 op¬ 
erators. All matter fields belong to lOj and 5j representations (leaving aside right-handed 
neutrinos in the singlet representation that may also be present), where i = 1,2,3 is the 
generation index. 

In order to connect the minimal renormalizable supersymmetric SU(5) model with ex¬ 
perimental data, one has to deal with three different theories: SU(5) above the unification 
scale Mguti the Minimal Supersymmetric Standard Model (MSSM) between Mgut and 
the supersymmetry scale nisusyi and the Standard Model (SM) between niausy and the 
weak scale rriz- Since the heavy GUT states (resp. the superpartners) are not degenerate 
in mass, the matching between the SU(5) theory and the MSSM at Mgut (resp. between 
the MSSM and the SM at msusy) will involve threshold corrections. 

In the next subsections we give the relevant parts of the corresponding Lagrangians 
(i.e. the Higgs and Yukawa sectors and the soft supersymmetry breaking terms, which 
determine the superpartner spectrum) and we specify our notations and assumptions. 

2.1 The SU(5) model 

The superpotential of the minimal renormalizable supersymmetric SU(5) model is deter¬ 
mined by its held content, gauge invariance and renormalizability. It can be divided into 
two parts describing the Higgs and Yukawa sectors, respectively: 

Wh = ^ Tr 24|^ -I- Tr 24f^ -I- 5//(/X// -I- 24//) 5//, (2.1) 

H^y = ^A^°10H0,-5//-v/2ApH0,-5//, (2.2) 

in which we have omitted terms involving right-handed neutrinos as well as R-parity vi¬ 
olating couplings that may be present, depending on how neutrino masses are generated. 
After having solved the equations of motion for the 24// in the SM singlet direction: 

(24//) = (To Diag(2,2, 2, —3, —3) , Uq = /xs/As, (2.3) 

and performed the fine-tuning needed to achieve doublet-triplet splitting in the Higgs sector, 
one can write down the masses of the heavy states in terms of the SU(5) superpotential 
parameters: 


mv = bgcurCTo , rriT = briH(^o, mg = -mg = -bmg = bX^a^ , (2.4) 

where my is the mass of the SU(5) gauge bosons in the representations (3, 2)_5/3©(3, 2 )+ 5 y 3 
of the SM gauge group; is the mass of the colour triplet and antitriplet pair (T, T) 
contained in 5// © 5//; and mi, mg and mg are the masses of the SM singlet, SU(2) triplet 
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and SU(3) octet components of 24/^, respectively. Demanding that the snperpotential 
conplings A 5 and rju be in the pertnrbative regime and taking into account the fact that 
the unihed gauge coupling Qquj. is of order 1 , one obtains the following constraint: 

rriT, ms, m 3 < my ■ (2.5) 

We also assume that supersymmetry breaking is coming from above the GUT scale, as 
for example in supergravity. In practice this means that the soft terms should be SU(5) 
symmetric at the GUT scale: 

- Csoft = {ml^)ij lolfo,- + {ml)ij + Q A]^ 10 . 10 ,-5^ - 4 §U0,-5^ + h.c.) 

+ {b^ 5h5h + h.c.^ + — Mi /2 A“A“. (2.6) 

We will consider the possibility of generation-dependent soft terms (as explained in the 
introduction, generation-dependent A-terms are needed to correct the SU(5) fermion mass 
relations), but in order to comply with the strong constraints coming from flavour physics 
we must ensure that they do not induce large flavour-violating effects. To this end, we 
assume that the soft sfermion mass matrices are diagonal in the basis in which the Yukawa 
couplings A® are diagonal, and that the A-term matrices are diagonal in the corresponding 
fermion mass eigenstate basis: 

("^10)0 = > 4 = 

AU = Al^5ij in the basis Ah’ = Al^5ij , ( 2 . 8 ) 

so that all flavour violation at the GUT scale is concentrated in the up squark sector and 
controlled by the GKM angles, yielding an effective alignment of sfermion soft terms with 
fermion masses [18]. In addition, we assume that the soft masses of the hrst two generations 
of sfermions are degenerate: 

^loi = ^lo2 > ■ ( 2 - 9 ) 

Finally, we will take M 1/2 and the A-terms (as well as the /i parameter /r = — Si/Hao) 

to be real. This may be more than what we need to evade flavour and GP constraints from 
low-energy experiments, especially in view of the fact that the superpartner spectrum is 
heavy, but this choice also helps reducing the number of parameters. In our subsequent 
exploration of the parameter space of the minimal renormalizable supersymmetric SU(5) 
model we shall completely neglect flavour violation in the sfermion sector, including the 
small amount of flavour violation that is generated from the running of the soft terms. 

2.2 MSSM 

Below the scale Mqut, the relevant theory is the MSSM, with snperpotential 

Wmssm = + fiHuHd, ( 2 . 10 ) 

and soft supersymmetry breaking terms 

- Esoft = {rnA)ij u-uf + 

+ Elqjhu - Af^ diqjhd - Afj e'^l/hd -h h.c.) 

+ + m]j^h\hd + {B Khd + h.c.) 

+ + ^M 2 WW + ^M3rr, (2A1) 
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where the contraction of SU{2)l indices is understood (for instance, HuHd stands for 
= HiHj - HlHl where i,j = 1,2 are SU (2 )l indices and etj is the totally 
antisymmetric tensor with eu = +!)• Due to the boundary conditions (2.7) and (2.8), 
and to the fact that we are neglecting the effects of the CKM matrix in the running, the 
sfermion soft terms keep a diagonal form all the way down to low energies. 

2.3 Standard Model 

Below the matching scale 

^^susy — \J ij^susy') i.^^susy') — \J ij^susy') {j^susy) (^D2) 

(where the last approximation is valid as long as the mixing in the stop sector is small, 
i.e. rritXt -C rriQ^^mlc), the relevant theory is the Standard Model. The Higgs potential 

is given by (with the SM Higgs doublet given by h = a/2 (cos (3 hd + sin {3 icr'^hl) in the 
decoupling limit): 

VsM =-mlh^h+, (2.13) 

while the Yukawa Lagrangian is 

^Yukawa ^ij Qj ^ij ^Ri^^Qj ^ij T h.C. , (2.14) 


where h = ia‘^h*. 


3 Renormalization group equations 

In this section, we collect the SM and MSSM renormalization group equations (RGBs) and 
various expressions used in our analysis (from boundary to matching conditions). We use 
the SM RGBs [28] between mz and rrisusy) and the MSSM RGBs [29] between rrisusy and 
Mgut- The gauge and Yukawa couplings, as well as the Higgs quartic coupling are evolved 
with the 2-loop RGBs, with the 1-loop threshold corrections accounting for the splitting 
of superpartner masses added at the scale rrisusy All soft parameters (A-terms, gaugino 
masses and soft scalar masses) are run at 1 loop. 


3.1 Gauge couplings 


The 2-loop RGBs for the gauge couplings gt {i = 1, 2, 3 for the gauge groups U{1)y, SU{2 )l 
and SU{3)c, respectively, with gi = ^5/3 g') read: 



(dvr)^ * (dvr)^ 



- 5] G„Tr(AtA„) 




(3.1) 







where t = ln(m/mo), m being the renormalization group scale, and the /3-function coeffi¬ 
cients below and above nisusy are given by: 


6f^ = (41/10,-19/6,-7), 


bfSSM ^ (33/5^ 


■ 199/50 27/10 44/5 ' 


■ 199/25 27/5 88/5 

tdSM _ 
^ij 

9/10 35/6 12 

tjMSSM _ 

9/5 25 24 


11/10 9/2 - 

-26 


11/5 9 14 


■ 17/10 1/2 3/2 ■ 



■ 26/5 14/5 18/5 ' 

r^SM _ 
^ia. 

3/2 3/2 1/2 

5 

r'MSSM _ 
^ia 

6 6 2 


2 2 0 



4 4 0 


(3.2) 

(3.3) 


(3.4) 


In the last term of Eq. (3.1), the MSSM Yukawa couplings should be replaced with the 
SM ones {Ha) below rrisusy 

At rrisusy, the running gauge couplings should be converted from the MS scheme to the 
DR scheme, in which the MSSM RGEs are written: 



[^susy 

= a-^ 

\^susy 

-e) 

5 


[^susy “1“ 

= a2^l 

\^susy 

-e) 

1 

Ovr 


[^susy 

= “3^1 

\^susy 

-e) 

1 

47r 


where ctj = gll{Hi) and e —)■ 0+. 


(3.5) 

(3.6) 

(3.7) 


3.1.1 Threshold corrections to gauge couplings 

Imposing gauge coupling unihcation at the GUT scale: 

Oil{MGUT) = 012{MgUt) = 0'3{Mgut) = OiGUT (3.8) 


implies certain relations among the masses of the various thresholds (supersymmetric part¬ 
ners of the SM helds and heavy GUT fields). Adding 1-loop threshold corrections [30, 31, 
32, 33, 34] to the running gauge couplings evolved with the 2-loop MSSM RGEs between 
the scales rususy and Mgut yields the following relations {i = 1, 2, 3): 


2-loop single-scale MSSM 


1-loop low-energy threshold corrections 


a 


-1 

GUT 


= (Mgut) — 


27r V 


^,(n) 


TTlr. 


1 

Hr L 


A6; In 


Mr- 


mv 


+ Abf In 


Mr 


1718 


+ In 


Mr 


1713 


+ Abf In 


Mr 


rriT 


l-loop high-energy threshold corrections 


(3.9) 


where the denote the values of the gauge couplings obtained by solving numerically the 
2-loop MSSM RGEs with all superpartner masses at the scale rrisusy, and n runs over the 
superpartners. Their contributions Ab'f^'^ to the /3-function coefficients are given by: 


Afef' 

= (1/30,1/2,1/3), 

A6? = 

(2/5, 2/3,0), 

Abf 

= (4/15,0,1/6), 

Abf = 

(1/10,1/6,0), 

Abf 

= (1/15,0,1/6), 

Abi = 

(0,0,0), 

Abf' 

= (1/10,1/6,0), 

Abf = 

(0,4/3,0), 

A6f 

= (1/5,0,0), 

A6f = 

(0,0,2), 
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(3.11) 


while the contributions of the heavy GUT fields are: 

Abf = (0,0,3), Abf = (0,2,0), 

Ab] = (0,0,0), Abf = (2/5,0,1), 

AbY = (-10,-6,-4). 

Taking appropriate combinations of the three equations (3.9), one obtains: 


rriT 

Mgut 


= exp 


X 


m 


3 1/3 


M, 


GUT 
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^GUT 



(3.12) 


(3.13) 


2ir 
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1/5 4/51 
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27r 


In 


19 T9 T 


m 


91 

susy 


n 

i=l 


IT T 

m^nrif 
m\cmlc 


1 1/18 


(3.14) 


At the 1-loop level, the matching scales Mgut and rrisusy drop out from Eqs. (3.12) 
and (3.13), while only rususy drops out from Eq. (3.14). Since rususy one can 

neglect the mixing between the higgsinos and the electroweak gauginos and identify: 

m-h = fi, (3.15) 

tjia = 


p mH^{msusy))il + 1/ tan^^) 

- \l{'mH^imsusy) - mH^{msusy)){taii^ -1-1) / (tan^/3 - 1) , (3.16) 

rriuj/mg ~ a2{mg)/as^rrig) , (3.17) 

where tan/? = (h.°)/(h°) is the ratio of the two MSSM Higgs doublet vevs, and p satisfies 
the electroweak symmetry breaking (EWSB) condition (again neglecting m|): 

tan^ /? 


|^|2 ^ ^HaK^susy, 


m 




m 


susy ) 


tan^ /? — 1 


(3.18) 


In spite of the initial historical success [35, 36, 37, 38], it is well known that gauge 
couplings do not unify accurately at the 2-loop level in the MSSM with TeV-scale super¬ 
partners. High-energy threshold corrections thus play a crucial role in achieving precise uni¬ 
fication [33]. Using Eqs. (3.12) and (3.13), one can express the combinations of GUT state 
masses needed for exact 2-loop unification in terms of the superpartner masses [33]. As¬ 
suming that all superpartners have masses equal to mgusy, one obtains for the colour triplet 
mass and for the combination of heavy gauge boson and adjoint Higgs masses m^^mg'^^: 

\ 5/6 

t \ 


rriT 


2 X 10^^ GeV 


2/3 


rriy m 


1/3 

3 


2 X 10^® GeV 


^^susy 

ITeV 

^^susy 

ITeV. 


) 

\ -2/9 


(3.19) 

(3.20) 
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where we have used the fact that m 3 = mg in the minimal renormalizable supersymmetric 
SU(5) model. For superpartner masses in the TeV range, the colour triplet is far too light 
and makes the proton decay too fast, which led Ref. [7] to conclude that the minimal 
renormalizable supersymmetric SU(5) model is excluded (this conclusion has been found 
to be mitigated at the three-loop level [39] though, and can be avoided for a specihc choice 
of the soft terms [3]). 

3.2 Yukawa couplings 

For the Yukawa couplings, we use the 2-loop Standard Model RGBs [28] below the scale 
rUsusy, and the 2-loop MSSM RGBs [29] above it. We neglect all GKM contributions and 
work with diagonal Yukawa matrices. In the absence of threshold corrections, the matching 
conditions at rrisusy are: 



/ sin, 

(3.21) 


COS, 

(3.22) 


hr^^^eij^susy) / COS/3 . 

(3.23) 


where the h couplings (resp. the A couplings) are the diagonal entries of the SM Yukawa 
matrices i/Fn'Y (^^esp. of the MSSM Yukawa matrices 


3.2.1 Threshold corrections to Yukawa couplings 

At the GUT scale, the SU(5)-invariant boundary conditions apply: 

^u{Mgut) = ^u{Mgut) , (3.24) 

^d{.Mgut) = ^'e{Mgut) ■ (3.25) 

Threshold corrections due to the splitting of the heavy GUT state masses slightly modify 
these relations (see later). After running down to low energy, Bqs. (3.25) lead to predictions 
for down-type quark and charged lepton masses that are in gross contradiction with the 
data. Supersymmetric threshold corrections at the scale mgusy may cure this problem. 


Supersymmetric threshold corrections to light fermion masses 

In the following, we will neglect supersymmetric threshold corrections to the leptonic 
Yukawa couplings, and consider only the corrections to the down-type quark Yukawa 
couplings®, whose dominant contributions are proportional to 03 and Xf (see however 
Ref. [40]). This will allow us to derive the SU(5) Yukawa couplings A® by simply running 
the charged lepton couplings up to the GUT scale. The leading supersymmetric threshold 
corrections to down-type quark masses are given by the gluino and higgsino contributions 
[ 8 , 9, 10, 12, 13, 14] (the latter can safely be neglected for the hrst two generations): 


Arud, 


20:3 

dvr 


m~gXd^ 




dR, 


lOvr^ 


nXt tan^/ 3 (|/ip,m? ,m|)5i3, (3.26) 


^Supersymmetric threshold corrections to up-type quark masses remain under control, as no large A- 
terms are needed to correct the SU(5) prediction. As for the top quark mass, even the large stop mixing that 
may be needed to reproduce the measured Higgs boson mass does not induce sizable threshold corrections. 
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where 


Xt = At !At - /i cot /? , 
Xd, = AdJXd, - /itan/3, 
Xei = Ale^Ae, -/itan/3, 

and the loop function J 3 is dehned by: 


(3.27) 

(3.28) 

(3.29) 


h{x,y,z) = 


with the limits 


h{x,x,z) = 


xyln {x/y) + yzhi {y/z) + zx\n. {z/x) 
(x - y)(y - z)(z - x) 

1 — (z/x) + (z/x) In (z/x) 


x(l — (z/x))'" 


l 3 (x,x,x) = — . 

lx 

The matching is done at the scale rususy'- 


(3.30) 


(3.31) 

(3.32) 


where 


m 


SM 

di 


™MSSM 



^d, ) 




= hd^(msusy)v, = Xdii'msusy) V cos(3 , 


MSSM 


(3.33) 

(3.34) 


in which v = (h.°). As a hrst approximation, b-r Yukawa unihcation is a relatively suc¬ 
cessful prediction of SU(5), while the discrepancy between the prediction and the data is 
much more important for the hrst two generations. As can be seen from Eq. (3.26), the 
non-holomorphic (cx /itan /3) contributions to Arud/md and ArUs/ms are the same for equal 
hrst two generation squark masses, while the ratios irts/m^ and md/me are widely diherent. 
This implies that large A-terms Ad and Ag are needed to bring these ratios into agreement 
with experimental data®, which in turn makes the electroweak vacuum metastable [12, 13]. 
This issue will be discussed in Section 6 . 


High-scale thresholds corrections to A;, and A,^ 

In addition to supersymmetric corrections at the superpartner mass scale, Yukawa cou¬ 
plings are also subject to high-scale threshold corrections due to the heavy GUT states. 
These may ahect in particular bottom-tan Yukawa unihcation, which as explained before is 
an important constraint on the model, hence one must take them into account. In practice, 
all one needs is the diherence between Xb(MGUT) and Xt-(Mgut) induced by the GUT-scale 
threshold corrections. One can check that it is given by: 


XhiMcur) — XriMcux) 


Xt(Mgut) 

(47r)2 


A?(M< 


GUT 


MguT . 2 1 AIguT 

In-In- 


rriT 


mv J 


(3.35) 


GUT threshold corrections also ahect strange quark-mu and down quark-electron Yukawa 
unihcation, but the numerical ehect is negligible compared with the size of the low-scale 
(supersymmetric) threshold corrections that are needed to account for the observed masses. 

^Incidentally, it turns out that the correct ratio cannot be obtained from the corrections pro¬ 

portional to /rtan/3 and to Xt alone, and that a large At, is also needed. 
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3.3 Higgs quartic coupling 


For heavy stop masses, the proper way to compute the lightest supersymmetric Higgs 
boson mass (for the standard computation, see Refs. [41, 42]) is to consider the effective 
theory below the scale mgusy, in which all superpartners and heavy Higgs bosons have been 
integrated out and the Higgs boson mass is determined from the Higgs quartic coupling A, 
with the value of X{mgusy) determined by the supersymmetric theory valid above msusy At 
tree level, the matching condition is A = {?)g\/b + 5 ^ 2 ) cos^ 2/3/4, while at the 1-loop level it 
is given by: 

X{mgusy) = i^^glimsusy) + glimsusy)) + AA^^" + AA^ + AA^ (3.36) 


where the hrst term on the right-hand side of Eq. (3.36) is the tree-level contribution, 
accounts for the conversion of the gauge couplings from the DR scheme to the 
MS scheme, and AA*^ and AA^ are the one-loop threshold corrections due to scalars and 
electroweak gauginos/higgsinos, respectively, whose expressions can be found in Ref. [43]. 
The dominant contributions are the ones proportional to (mt/n)^ in AA*^, which in the case 
where all sparticle masses lie close to rrisusy (and in particular ttiq^ ~ ~ msusy) reduce 

to the leading stop mixing term: 


AA'^ ~ 


(47r)2 
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(3.37) 


For large values of tan/3 and/or large values of one should also include on the RHS 
of Eq. (3.37) the leading sbottom and stau contributions, which in the limit mg^ ~ ~ 

~ mgg ~ msusy read: 

Qht XI r 1 XI 

{\'Kymimi^ [ Vlmi^mi^ 


+ 


2ht 


XI 


(dvr)2 mr^mf^ 


1 - 




12 mf^mf^ 


(3.38) 


Considering only the leading term (3.37), one can see that for each value of msusy 
and tan (3 there exist either 4, 2 or 0 different values of Xt satisfying the matching condi¬ 
tion (3.36). There also exists a lower bound on msusy for each value of tan /3, reached when 
Xt becomes 

Xr" = ±V6 msusy, (3.39) 

where the stop mixing contribution reaches its maximum [44] (when higher order corrections 
are included the Higgs mass also depends on the sign of Xt). This is illustrated in Fig. 2, in 
which the lower bound on mgusy is represented by the green curve, assuming a superpartner 
spectrum as in Fig. 1. The comparison with the red curve shows the importance of the 
1-loop threshold corrections to the Higgs quartic coupling. Note that the SM parameters 
were set to their central values in this hgure; taking into account the uncertainty on the top 
quark mass would spread the curves into bands which become very broad at small tan/3 
values. For a discussion on this point, see Ref. [44]. 


3.4 A-terms 

As discussed in the introduction and in Subsection 3.2, we assume the A-term matrices to 
be diagonah in the corresponding fermion mass eigenstate basis, with generation-dependent 

^Due to the CKM mixing in the quark sector, this assumption is not renormalization group-invariant. 
However, the off-diagonal entries generated by the running from the GUT scale to low energy are suppressed 
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tanjS 


Figure 2: Higgs mass constraint in the (tan /?, rrisusy) plane for different matching conditions, 
with all SM parameters set to their central values. The red curve corresponds to the tree- 
level matching condition and the blue curve to the full 1-loop matching condition (3.36), 
while for the green curve only the leading 1-loop threshold correction (3.37) was used. 
Both the blue and green curves assume maximal positive stop mixing (X* = VGrUgusy)- 
The superpartner spectrum is chosen as in Fig. 1. 

entries: 

He/ = Diag(H„,Hc, A), (3.40) 

Ad = Diag {Ad, H*, Ab ), (3.41) 

Ae = Diag (He, H^, H^). (3.42) 

These matrices are run with the 1-loop MSSM RGBs, with the SU(5) boundary conditions 
imposed at the GUT scale: 

Au{Mgut) = Ae{Mgut) , (3.43) 

Ad{Mgut) = A'^{Mgut) ■ (3.44) 

3.5 Gaugino masses 

For the running of gaugino masses we use the 1-loop RGBs: 

= (3.45) 

at zn 

and impose the SU(5) boundary condition at the GUT scale: 

Mi{Mgut) = M2{Mgut) = ^^{Mgut) = ^ 1/2 • (3.46) 

by the small CKM angles (the same statement holds for the soft sfermion mass matrices). We shall therefore 
neglect this effect in the following. 
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3.6 Soft scalar masses 

We also assume the soft sfermion mass matrices (X = L,d^) to be diagonal 

with generation-dependent entries: 

m\ = Diag(m|^,m^^,mU , (3.47) 

however with imposed at the GUT scale. Furthermore, SU(5) invariance 

requires the following relations to hold at Mgut ii = 1,2,3): 

'i^qS^gut) = ^uiiMcur) = '^e<r{MGUT) = , (3.48) 

’’T^Lii^Gur) = mj^c^McuT) = rn^i- (3.49) 

Hence the splitting of the soft sfermion masses within SU(5) representations and between 
the hrst two generations is only due to the running, performed at the 1 -loop level as for the 
other soft terms. As for the soft Higgs masses, we allow the possibility of different boundary 
conditions for the two Higgs doublets of the MSSM, namely rn?^^{MQUT) 7 ^ rri^j^^McuT)- 
Note that our boundary conditions are less restrictive than the so-called minimal super¬ 
gravity ansatz (mSUGRA), which assumes universal scalar and gaugino masses as well as 
A-terms proportional to the Yukawa couplings. By contrast we allow for some generation 
dependence in the sfermion soft terms, but we require them to be aligned with fermion 
masses (exactly at the GUT scale, approximately at low energy due to the running) in 
order to minimize supersymmetric contributions to flavour-violating processes. 

3.7 /i and B^i terms 

Gontrary to the other MSSM parameters, the /i and Bjji terms are not hxed at the GUT scale 
and renormalized down to low energy, but rather determined from the tree-level electroweak 
symmetry breaking conditions at the scale mgusy, as is usually done (by inserting the 
running soft terms in the tree-level Higgs potential one ensures that the most relevant 
1-loop radiative corrections are taken into account). This procedure yields: 

I |2 _ rn^H^ - f3 m| 

tan2/?-l 2 ’ 

^ / 9 9 , , 9 x sin 2/3 

B = (m^^+ 2|/i| ) 

-mjjJtan/3 _ m|tan/3 
tan^ /3 — 1 tan^ /3 -f 1 ’ 

while the sign of fi remains undetermined. 

4 Solutions of the RGEs for soft terms 

Here we present (approximate) semi-analytic solutions to the 1-loop renormalization group 
equations for the MSSM soft terms. More details about the procedure used to derive them 
can be found in Appendix A. 


(3.50) 

(3.51) 
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4.1 Gauge and Yukawa couplings 

As a first step, for each choice of rrisusy and tan/3, one runs the gauge couplings gi(t) {i = 
1, 2, 3) and the Yukawa couplings A„(f) (n = t, b, r, c, s, /i, u, d, e) from tz = Inm^ to ts = 
InrUsusy with the 2-loop SM RGBs. Then, after applying the matching conditions (3.5)- 
(3.7), (3.21) and (3.23) at the scale rngusy) fhe gauge couplings and the up-type quark and 
charged lepton Yukawa couplings are further evolved up to the GUT scale with the 2-loop 
MSSM RGBs. At Mcut the SU(5) boundary conditions (3.25) and the heavy threshold 
corrections (3.35) are imposed, then the down-type quark Yukawa couplings are run back 
from fo = InMct/'T to ts = \^msusy (in this procedure one needs to hx the heavy GUT state 
masses my and mj^, in addition to the GUT scale itself). This provides us with numerical 
solutions for the running gauge couplings gi{t) {i = 1,2,3) and Yukawa couplings A„(f) 
{n = t, b, T, c, s, /i, M, d, e) in the range ts <t < to, accurate at the 2-loop level. 


4.2 Gaugino masses 

The solutions to the 1-loop RGBs for gaugino masses read: 

M,(m) = My, e- . (4.1) 

These 3 equations connect linearly 4 variables: Mi^ 2 , 3 {,fnsusy) and Mi/,. 


4.3 A-terms 


As explained in Appendix A, the hierarchy of Yukawa couplings makes it possible to solve 
the 1-loop RGBs for A-terms in a sequential way, from At to Ag. With the convention that 
the index n runs over the ordered values {t, b, r, c, s, /i, u, d, e} (such that e.g. p < c means 
p = t,b or t), one can write the (approximate) solutions as: 


An{t) 


A„(4)d‘. I''' 


+ 






(4.2) 


where is a function of the gauge and Yukawa couplings, and y/f is the product of A„ 
times a linear combination of gfMt and XmAm with m < n, all of which are already known 
quantities, since the RGBs for the A^’s are solved in order of increasing m. One can 
therefore rewrite Bq. (4.2) in the form: 


Mt) = E ^nm {t)Amito) + bn{t)Mi/, , (4.3) 

m=l 


in which the coefficients anm{t) and bn{t) are integrals that can be evaluated numerically 
after having solved the MSSM RGBs for the gauge and Yukawa couplings. Imposing the 
SU(5) boundary conditions Ad-(to) = Ae^{to), one obtains the running A-terms at an arbi¬ 
trary scale as linear combinations of the SU(5) soft parameters Mi/,, A]^, A\ {i = 1,2,3), 
with numerical coefficients depending on the choice of mgusy and tan (d (and very mildly on 
Mqut, 'rriT and my). 

In practice it may be more convenient to express the An{t) in terms of other input 
parameters, for example Moijrisusy) and the AdX'^susy) {i = 1,2,3), since these quantities 
enter the supersymmetric threshold corrections to the down-type quark masses needed to 
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fit the experimental values. One may also want to trade for At{msusy), which is directly 
constrained by the measured value of the Higgs mass. To this end, it suffices to invert the 
linear relations (4.3) so that the can be rewritten as a function of the desired input 

parameters. In the following, we shall choose: 

Al’^irfigusy) 1 An^irfigy^gy^^ Aj^^l^uigy^gy') (i 1,2,3). (4.4) 

4.4 Soft scalar masses 

Following Appendix A, we hrst introduce the combinations of masses (which appear in the 
RGBs for soft scalar masses): 


3 


s = ^ [m|^ - - 2m| + m| + m|] , 

i=l 

(4.5) 

and 


2 1 2 1 2 
= rriQ^ + m^c + tuh^ , 

(4.6) 

2 1 2 1 2 

Sd, = +m^^. 

(4.7) 

= ml, + m% + . 

(4.8) 

The 1-loop RGB for S is easily integrated to give: 


s(t) = 

(4.9) 


where S'(to) = '^r/„(0) — ''^^■^(0) due to the SU(5) boundary conditions on soft scalar 
masses, while the RGBs for the S„’s can be solved sequentially in a similar way to the A- 
term RGBs, taking advantage of the hierarchy of Yukawa couplings to neglect subdominant 
terms. We then dehne: 

= 7^ fdt'Alit'), (4.10) 

(47r)"= Jto 

^ ttY /' dt'gUt')Mm , (4.11) 

(47r)"= Jto 

Is{t} ^ dt'gUt'}S(t'). (4.12) 

(47r)^ Jto 

With these ingredients one can express the solutions to the 1-loop RGBs for soft scalar 
masses as: 

<(t) = rnKu) + (E„(«) - E,.(io)) + E + E KImM + R^hW , (4.13) 

n n i 

where the index a runs over Qsgp, 2 u ^3 2 d-^ 3,2,i) ^^3 2 i}) numerical 

coefficients T“, f/" and 14“ (resp. i?“) are given in Table 4 (resp. Table 3) of Appendix 
A. Plugging the previously derived semi-analytic expressions for A„(t), Bq. (4.3), and for 
E„(t) into Bq. (4.13), one then arrives at 

+ Y.danm{t)An{to)Am{to) 

P n,m 

+ X/ ^ctn{t)An{to)M i/2 + fa{t)Mi/2 , (4-14) 
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in which the coefficients Ca/sit), danmit), Canit) and fait) are integrals that can be evaluated 
numerically just from the knowledge of the solutions to the 2-loop MSSM RGBs for gauge 
and Yukawa couplings. Finally, imposing the SU(5) boundary conditions on soft terms, one 
obtains the running soft scalar masses at an arbitrary scale as quadratic functions 

of the SU(5) soft parameters miQ., mg., mn^^McuT), ^ 1 / 2 , and {i = 

1,2,3), with numerical coefficients depending on the choice of nisusy and tan/? (and very 
mildly on Mguti rriT and mv)- 

As was done for A-terms, one can easily rewrite the running soft scalar masses in terms of 
other input variables by inverting the system of equations (4.3) and/or (4.14). As explained 
before, a convenient choice for exploring the parameter space of the minimal renormalizable 
supersymmetric SU(5) model is to trade the 7 SU(5) parameters M 1 / 2 , and A/ for 
Ms^nisusy), Adi{msusy) and Aui(msusy)- One may also invert 8 of the 17 equations (4.14) 
in order to replace the GUT-scale masses mio^, mg., mH^^Mcur) and mn^iMcuT) by 8 
low-energy soft scalar masses, so that (for given values of rngusy and tan/3) the whole 
supersymmetric spectrum is parametrized by 15 low-energy input variables. 


5 Proton decay 

Proton decay is one of the main prediction of Grand Unihed Theories, and since it has 
not been observed yet, it sets strong constraints on the parameter space of the minimal 
renormalizable supersymmetric SU(5) model. Qualitatively, the proton lifetime behaves as: 


Tp ~ t{p ^ K^v) oc mym^„gj^tan^/3/(l + tan^/3)^ , (5.1) 

implying a tan /3-dependent lower bound on the superpartner mass scale. 

To compute precisely the proton lifetime we will need the following input parameters: 


mp = 0.9383 GeV, 
03(2 GeV) = 0.31, 
GeV) = (0.111 ± 0.027) GeV^, 


= 0.4937 GeV, (5.2) 

a 3 (mQ = 0.22, (5.3) 

lVo^2^(2GeV) = (0.036 ± 0.014) GeV^ , (5.4) 


where and appear in the hadronic matrix elements for proton decay, as well as 
the entries of the GKM matrix: 


VcKMijnz) 


1 — A^/2 A AX^{p — iri) 

-A 1 - AV2 AA2 

AA^(1 —p —ir;) —AA^ 1 


(5.5) 


here written in the Wolfenstein’s parametrization, with [45]: 


A = 0.22537, A = 0.814, p = 0.117, p = 0.353 . (5.6) 


Details about the computation of the proton lifetime can be found in Appendix B. The 
predicted proton lifetime is compared with the experimental constraint r(p —)■ K^v) > 
2.3 X 1033 yrs (90% G.L.) [46]. 


6 Vacuum (meta) stability 

In the general MSSM, some regions of the parameter space lead to instabilities of the 
electroweak vacuum. One may encounter two kinds of dangerous situations. 
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The first one is the possible existence of directions in held space along which the potential 
is unbounded from below (UFB). To remain on the safe side, we will allow only the points 
in parameter space that do not possess any such direction. The associated constraints on 
the model parameters will be summarized in Subsection 6.1. 

The second one is the metastability of the electroweak vacuum, which due to the large 
trilinear soft terms needed to correct the SU(5) predictions for fermion masses cannot be 
avoided. This means that there are minima in held space, lower than the electroweak vac¬ 
uum, into which it will eventually decay. These new minima are the so-called charge and 
colour breaking (CCB) vacua. Although one cannot forbid the decay of the electroweak 
vacuum into a lower minimum, one can check whether its lifetime is long enough on cosmo¬ 
logical time scales. The procedure followed in our analysis is summarized in Subsection 6.2. 

All computations in this section are done at the tree level. The generalization to higher 
orders is conceptually straightforward, although technically much more involved, so we will 
leave it for future work. Below we summarize the discussion presented in Appendices C 
and D, where technical details and relevant references can be found. 

6.1 Unbounded from below directions 

As explained in Appendix C, the tree-level constraints associated with the absence of UFB 
directions can be written as (neglecting m^): 

( 6 - 1 ) 

which must be satished for any of the three slepton generations. 

6.2 Charge and colour breaking vacua 

We shall discuss the constraint applying to At separately from the ones applying to all other 
A-terms. The reason to treat them differently is that in the second case the D-terms can 
be considered to be vanishing to a good approximation, thus providing constraints on the 
helds, while in the hrst case this assumption is not justihed due to the large top Yukawa 
coupling®. We will nevertheless set the colour D-terms to zero for simplicity, but allow for 
non-vanishing hypercharge and SU(2) D-terms. 

6.2.1 Constraints on A„ {n ^ t) 

Let us hrst dehne: 

M2(^, H, H) ^ ^ ^2^2 ^ ^ |^|2^ ^2 ^ ^ |^|2^ ^2 

H , H , c ,) = 2 ( 2^/>2 + ^ ^2 ^ ^ 2 ^ ^ 

and write the SU(2) D-term constraint: 

0 = + (3 - 2c)i)'^ + H^- H^, (6.5) 

®Strictly speaking, the same comment applies to At, in the large tan/? regime. In the minimal renor- 
malizable supersymmetric SU(5) model, however, large values of tan/3 are excluded by a combination of 
constraints (see Section 7), so only the case of a large At needs to be discussed separately. 


-2BHH, (6.2) 

(6.3) 

(6.4) 


19 



where "0, 'll! are sfermion fields and if, if Higgs fields that parametrize a specific direc¬ 
tion in field space along which a CCB minimnm is present. The different possibilities for 
the constant fields if}, ip, H, if, their mass parameters and the coefficients c and are 
summarized in Table 5 of Appendix D. 

To evaluate the lifetime of the electroweak vacuum, we consider the (normalized) bounce 
action, which for An 7 ^ At can be approximated by: 

5 = ^(^,^,if,F). (6.6) 

This action is then minimized by varying the direction in field space (V’, ip, if, if), subject 
to the constraint (6.5). In order for the lifetime of the electroweak vacuum to be larger 
than the age of the universe, the minimum of S must satisfy 

> 9 . (6.7) 

Notice that one needs to minimize S with respect to two variables only. Indeed, one of 
the variables ip, ip, if, if is fixed by Eq. (6.5), and the quantity ( 6 . 6 ) only depends on ratios 
of fields. The minimization is performed numerically. A point of the parameter space is 
admitted if the bounce action satisfies Eq. (6.7). 


6.2.2 Constraint on At 

We now define: 

+ + + + ( 6 . 8 ) 

M 3 = —2t‘^(AtHu — Xt^Hd), (6-9) 

* = 2(2e + Hl + Hl), ( 6 . 10 ) 

A = + + (6.11) 

Due to the large value of the top quark Yukawa coupling, the bounce action can no longer 
be approximated by Eq. ( 6 . 6 ). One should instead minimize 


S = 1 + 


/(«)-/(0)\ z^M? 


m 


Ml ’ 


with 


/(«) = 


7r^/6 


16.5 


+ 


(1 — 4k)^ (1 — 1 — 4/t: ’ 


M^ 


( 6 . 12 ) 


(6.13) 


The minimization goes again over two variables (for example Hu/t and Hd/t), and is done 
numerically. 


7 Results and discussion 

We are now ready to address the question we asked in the introduction, namely whether 
the minimal renormalizable supersymmetric SU(5) model can be considered as a viable 
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extension of the Standard Model. To answer this question, we must scan over the parame¬ 
ter space of the model and search for points passing all phenomenological and theoretical 
constraints® (precise gauge coupling unihcation, correct predictions for the Higgs boson 
and charged fermion masses, proton lifetime, experimental lower bounds on superpart¬ 
ner masses, metastability of the vacuum and perturbativity of the model). This is not a 
straightforward task, as the model involves a large number of parameters with boundary 
conditions dehned at the GUT scale, while most constraints apply at low energy. In order 
to ease the whole analysis, we will use the semi-analytic solutions to the soft term RGBs 
obtained in Section 4. This will enable us to perform a much more efficient scan - even 
though it still involves a large number of parameters. 

Before going into this programme, let us hrst try to identify the region of the parameter 
space in which viable points are likely to be found. For a qualitative discussion of how each 
phenomenological or theoretical requirement constrains the model, we will consider only 
two parameters, tan/? and the overall superpartner mass scale rrtsusy'- 

• a powerful constraint on the parameter space comes from accommodating the mea¬ 
sured Higgs mass, which provides a tan /?-dependent lower bound on nisusy (the lower 
the value of tan/?, the higher the value of msusy)- Sizable mass splittings among 
superpartners can change this bound, but not very drastically. 

• another important constraint comes from the non-observation of proton decay, which 
provides another lower bound on rrisusy, with a different dependence on tan /?. The 
actual value of the bound depends on the details of the superpartner spectrum. 

• requiring the masses of the heavy GUT states derived from gauge coupling unihcation 
to remain below the cut-off scale Mpianck bounds rrisusy from above. The main con¬ 
straint comes from the Higgs colour triplet mediating proton decay, whose mass rriT 
strongly depends (through the gauge coupling unihcation condition) on the higgsino 
mass /i. 

• the parameter space is also bounded by vacuum metastability constraints associated 
with large values of the A-terms. These are unavoidably present due to the sizable 
threshold corrections needed in the down-type quark sector. Furthermore, a large At is 
necessary to accommodate the observed Higgs mass in some regions of the parameter 
space. Although the required values of the A-terms do not necessarily threaten the 
metastability of the electroweak vacuum, signihcant regions of the parameter space 
where the lifetime of the universe would be too short are excluded. 

• the requirement that the third generation Yukawa couplings should remain perturba¬ 
tive up to the GUT scale excludes some portions of the parameter space in the small 
tan /? region (top quark) and potentially also in the large tan /? region (bottom quark 
and tan Yukawa couplings). 

• the perturbativity of the parameters of the SU(5) superpotential, reflected in the 

condition mT,m3,ms can easily be satished by using the freedom allowed by 

gauge coupling unihcation. Indeed, for hxed Mqut and superpartner masses, rriT and 

®As discussed in the introduction, dark matter and neutrino masses can be accounted for by separate 
sectors (in particular when the lightest neutralino is not a suitable dark matter candidate), so we do not 
include them in the list of constraints to be imposed on the model. 
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Figure 3: A naive estimate of the allowed parameter space of the minimal renormalizable 
supersymmetric SU(5) model {pale green region) in the approximation described in the 
text, where all sfermions are assumed to have the same mass rususy 

mym^ are constant (using the fact that the colour octet and weak triplet components 
of the adjoint Higgs field have equal masses, m 3 = mg), see Eqs. (3.12) and (3.13). 
One can therefore increase my and decrease m 3 = mg by diminishing A 5 and /is in 
the superpotential ( 2 . 1 ) in such a way that uo = /is/As increases. Simultaneously, pn 
should be made smaller so that m^ = St/hUo stays constant. 

One can illustrate how these constraints restrict the parameter space of the model 
by plotting them in the (tan/?, rUsusy) plane, in the spirit of Fig. 1 in the introduction. 
To be able to do this, we assume a simplified superpartner spectrum with a common 
scale fUgusy for all sfermion masses, the /i parameter, and the SU(5) gaugino mass 
parameter Mi/ 2 - The down-type squark A-terms and At are chosen to fit the fermion 
masses and the Higgs mass, while the slepton A-terms are taken to be 2/3 of their down- 
type quark counterparts^^, and Au and Ac are assumed to vanish. With m^^ and the B 
parameter determined by the electroweak symmetry breaking conditions, the whole MSSM 
spectrum, including the heavy Higgs masses, can be computed. While these inputs are not 
fully consistent with SU(5)-invariant boundary conditions at the GUT scale, they make it 
possible to display all the constraints discussed above in a two-dimensional plot. The result 
can be seen in Fig. 3, a simplihed version of which was shown in the introduction. Let us 
now comment it: 

• the measured Higgs boson mass excludes the region below the red curve (which cor¬ 
responds to maximal stop mixing); 

• the region below the green line is ruled out by the experimental lower bound on 
the proton lifetime (assuming the high-energy phases appearing in the proton decay 
amplitude, see Eq. (B.7), all vanish). The shape of this line can be understood by 

^*^This empirical factor roughly mimics the effect of the running from Mgut, where the SU(5) relations 
Ad, = Aei hold, to rrisusy 
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noting that the proton lifetime approximately scales as 

tIp^K'^u) oc mymg„gj^tan^/3/(l + tan^/3)^ , (7.1) 

while given the assumptions made on the superpartner spectrum, the colour triplet 
mass is proportional to^^ 

rriT oc (tan^ /3 + , (7.2) 

giving 

t{p ^ K+i?) oc mini ^ ^ ; ('^•3) 

• the region above the purple curve is excluded by the vacuum metastability constraint 
on At, while the region above the brown curve (which is almost entirely due to Aj,) is 
ruled out by the constraints on all other A-terms. In drawing these curves we used 
the metastability conditions derived in Appendix D, but they can be approximated 
to a very good degree by: 

|7lt| < , \Ab\ < + m^ , (7.4) 

for At and Ab, respectively. Note that the constraint on \Ab\ is more easily satished 
for large tan (3 values due to m^j^ = ^lusyi‘^ ~ 1) j 

• perturbativity of the top quark Yukawa coupling excludes the region to the left of the 
blue curve (which corresponds to X^Mcut) = !)• There is no similar constraint for 
Aft and A^, which are always in the perturbative regime in the region of the parameter 
space shown in Fig. 3; 

• the constraint mr < Mpianck/^^ (resp. mp < Mpianck) rules out msusy values above 
the solid black line (resp. the dashed black line), while there are no additional con¬ 
straints from m^,m^,mT < my < Mpianck/^O] 

• the pale green area is the region of the parameter space that is allowed by all the 
above constraints. 

The tentative conclusion one can draw from Fig. 3, even though the assumptions made 
on the soft terms are not consistent with SU(5)-invariant boundary conditions at the GUT 
scale, is that viable points satisfying all phenomenological and theoretical constraints are 
likely to be found in the region bounded by 2 < tan /3 < 5 and 100 TeV < m^usy ^ 1000 TeV. 
This dehnitely needs conhrmation from a more careful investigation of the parameter space. 
Namely, one should scan over the parameters of the model, which besides tan/?, sign(/r) 
and the heavy state masses include 15 soft terms: 8 sfermion (ruioi j.si ’^ 5123 ) Higgs 
(ms, mg) soft masses, 1 gaugino mass parameter (M 1 / 2 ) and 6 A-terms (A^° 2 , 3 ) 

Since these soft parameters are dehned at the GUT scale, one needs to run them down to 
the scale msusy, where most constraints apply. In order to simplify the problem, we shall 
use the semi-analytic approximate solutions to the soft term RGBs derived in Section 4. 
This will allow us to trade the GUT-scale soft parameters for low-energy ones (namely 

^^Indeed, given the choice gX = one has from the electroweak symmetry breaking condi¬ 

tions = ms„j,y(2tan^/3 —1) and = 2ml,j^sy{tarA j3+l). Plugging this, together with the assumption 
of equal sfermion masses, into Eq. (3.12), one arrives at rriT oc msl)sy(tan^ /3 -|- 1 )^T 2 ^ 
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the nisusy values of the running parameters M 3 , Au^c,t, and of 8 suitably chosen soft 

scalar masses) and to use the constraints at rrisusy to effectively reduce the number of free 
parameters. 

Let us describe more precisely the procedure that we are going to employ: 

1 . hrst choose a random point in the (3, rrisusy) plane (where rrisusy is now the 
matching scale between the SM and the MSSM), together with sign(/i) and some 
sensible values of Mqut, fnr and mv; 


2. solve numerically the 2-loop MSSM RGBs for the gauge and Yukawa couplings (tak¬ 
ing into account the GUT-scale relations (3.25) and the heavy threshold corrections 
to Yukawa couplings (3.35), as explained in Subsection 4.1). Then, following the pro¬ 
cedure of Section 4, express the running soft terms at an arbitrary scale as algebraic 


m:- 


and 


combinations of the rrisusy values of M 3 , Au^c,t: A^^sfi^ ^ 
which through the EWSB condition (3.50) can be written as a linear combination of 


m% , 

^1,2,3 


m 


Hu 


{rrisusy) and mjj^{m susy\ 


3 . impose the following 8 constraints on the rrisusy values of the soft terms: (i) the deh- 
nition of the matching scale rrisusy = \JmQ^{msusy)mu^{rrisusy)', (ii)-(iii) the equal¬ 
ity of the hrst and second generation soft sfermion masses, which for simplicity 
we impose at the scale rrisusy rather than Mgut'- rn?Q^{msusy) = rn?Q^{msusy) and 
m\{msusy) = rn"^^{rrisusy)] (iv) the gauge coupling unihcation condition (3.12), in 

which mrp and Mgut are hxed and /r and rriA are functions of rri\j^ and through 
the EWSB conditions; (v) the 1-loop matching condition for the Higgs quartic cou¬ 
pling (3.36), keeping as a hrst approximation only the leading term (3.37); (vi)- 
(vinjihe supersymmetric threshold corrections (3.26) needed for the down-type quark 
masses to match their measured values. Whenever these constraints admit several 
solutions, we explore all of them. Imposing them allows us to express 8 of the 15 
input parameters at the scale rrisusy ia terms of the remaining 7 ones and of already 
known quantities; in practice we choose the remaining 7 free parameters (the ones to 
which we assign random values in order to explore the model parameter space) to be 
nr^{rrisusy)■} rn^c ^{rrisusy)^ -M^{rrisusy) nnd Au^c{msusy)- 


4. for each point of the parameter space, dehned by the chosen values of rriQ^ ^{rrisusy)) 
n^diJ<nTsusy), M^{msusy) and Au,c{msusy) (in addition to the values of tan^, rrisusy, 
Mgut, rriT and my, already hxed in the hrst step), one can improve the analysis by 
running the soft terms again, this time with the full RGEs, and use the complete 
expression for the Higgs mass. Once this is done, the phenomenological viability of 
this point must be further checked: proton lifetime, consistency of the superpartner 
mass spectrum with the experimental limits, absence of UFB directions and metasta- 
bility of the electroweak vacuum (which is checked following the procedure described 
in Appendix D), perturbativity constraints. 


The procedure described above was applied by testing each randomly chosen point in the 
{tdHi (3, rrisusy) plane for 1000 diherent random conhgurations of the superpartner spectrum, 
which were obtained by scanning over rnQ_^ ^{rrisusy), rn^c ^{rrisusy) and M 2 ,{Mgut) in the 
range [Q A rrisusy,‘^rrisusy], and over Au,c{msusy) in the range [-2 rrisusy,‘^rrisusy]- The results 
of this exploration of the parameter space of the minimal renormalizable supersymmetric 
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Figure 4: Left panel; successful {green dots) and unsuccessful {red stars) points among 
1000 randomly selected points in the {tan (3, uisusy) plane. Right panel: successful points 
only, where the colour of a given dot represents the number of superpartner spectrum 
conhgurations respecting SU(5) symmetry at Mgut that passed all phenomenological con¬ 
straints (the number of solutions found increases from light yellow to dark green). Also 
shown is a particular solution {magenta dot) characterized by a low value of the matching 
scale rususy = ^J'nrQ^{msusy) ^^^^{msusy) and described in greater detail in the text below. 
In both panels, the red, green, purple and black lines represent the naive constraints from 
Fig. 3 associated with the measured Higgs mass, the proton lifetime, the metastability of 
the electroweak vacuum and the condition < Mpianck/^0, respectively. 

SU(5) model can be seen in Fig. 4. The plot in the left panel displays all the points of the 
{tan/3,msusy) plane that have been checked; the ones for which at least one conhguration 
of the superpartner spectrum passed all constraints are in green, while the ones that failed 
this test are in red. The plot in the right panel shows how frequently a solution was found 
in the model parameter space for each green point. The more vivid the colour of the point, 
the more spectrum conhgurations survived all the constraints. It turns out that the naive 
estimate of the allowed parameter space (the green region in Fig. 3), despite not respecting 
the SU(5) symmetry, tells us something about how likely it is for a randomly chosen point 
in the parameter space with given values of tan f3 and rUsusy to be compatible with all 
experimental and theoretical constraints discussed at the beginning of this section. In the 
part of the green region of Fig. 3 around several hundreds of TeV (close to the upper limit 
given by the requirement mp < Mpianck/^^)-, a randomly chosen point in the parameter 
space has a few percent probability of success, while the probability decreases when one 
exits the green region. One can see that the naive constraints associated with the Higgs 
mass {red line) and with the proton lifetime {green line) are very robust, which in the 
case of the Higgs mass can be traced back to the fact that the leading term in the 1-loop 
matching condition for the Higgs quartic coupling, Eq. (3.37), only depends on niausy in 
the maximal stop mixing case. The naive vacuum metastability bound {purple line) is 
less robust but remains a reasonable approximation to the exact condition, contrary to the 
naive mp < Mpianck/^0 constraint {black line). The main reason for this is that the black 
line assumed p,{msusy) = 'msusy, while a sizable portion of the viable points of the parameter 
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space feature smaller values of /r, thus effectively decreasing the value of rriT by virtue of 
the gauge coupling unihcation condition (3.12). 

While these results conhrm the naive expectation from Fig. 3 that the allowed parameter 
space of the minimal renormalizable supersymmetric SU(5) model is restricted to the region 
of heavy superpartners, isolated points with lighter supersymmetric particles are likely to 
be missed in this random search. In fact we were able to hnd a viable point in the parameter 
space with nisusy = 63 TeV and some squarks as light as 22 TeV, depicted by a magenta dot 
in Fig. 4. We give below the input parameters at the GUT scale as well as the values of all 
soft terms at the scale nisusy'- 


INPUT: 


Point in the {ta,n (3, nisusy) plane: 



tan (3 = 

: 4.0, 

(7.5) 

‘^susy 

: 63.0TeV, 

(7.6) 

Mgut = 

: 2.0 X 10^®GeV, 

(7.7) 

sign(/r) = 

■ + 5 

(7.8) 

Soft terms at Mgut: 



Au(Mgut) = 

Diag (2 X 10■^ 9 x 10“^ 13.9) TeV , 

(7.9) 

X^di^GUr) = ^e(MGUT) = 

Diag (-0.4,41.2,121.7) TeV, 

(7.10) 

Mi/2 = 

34.4 TeV, 

(7.11) 

^Hu{XXIgut) = 

106.2 TeV, 

(7.12) 

^Hu{MgUt) = 

531.2 TeV, 

(7.13) 

mio, = 

Diag (79.9, 79.9,125.7) TeV, 

(7.14) 

= 

Diag (151.8,151.8,141.1) TeV. 

(7.15) 


OUTPUT: 

Values of gauge and Yukawa couplings at nigusy? in the supersymmetric 
{^susy + e) and in the non-supersymmetric theory (nisusy — ^)' 


9i(msusy + e) = (0.47777,0.61950,0.88966), (7.16) 

gi(msusy-€) = (0.47777,0.61900,0.88742), (7.17) 

{Xt,XbAr)imsusy + e) = (0.74778,0.06325,0.04229), (7.18) 

(ht,hb,K)(nisusy-e) = (0.72553,0.01123,0.01026), (7.19) 

(Ae, A„ Xy)(nisusy + e) = (2.61 x lO’^, 3.95 x 10-^ 2.52 x lO’^), (7.20) 

(K, hs, hf,){nisusy - e) = (2.53 x 10■^ 0.23 x 10■^ 0.61 x 10“^), (7.21) 

(A„, Ad, X,)(nisusy + e) = (5 x 10-^ 19 x lO"®, 12 x lO"®), (7.22) 

(hu, hd, he){msusy " e) = (5 X 10-^ 12 x 10“®, 3 x 10“®), (7.23) 
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Soft terms at m, 


susy 


-^u ij^susy) 

= Diag(0,0,-25.9) TeV, 

(7.24) 

-^dij^susy^ 

= Diag(-0.9, 87.6, 237.2) TeV, 

(7.25) 

-Ag ijTlsusy^ 

= Diag (—0.6, 55.8,163.4) TeV, 

(7.26) 

Mi{mg) 

= Diag(16.9, 29.9, 62.0) TeV, 

(7.27) 

ij^susy^ 

= 82.1 TeV, 

(7.28) 


= 483.5 TeV, 

(7.29) 

^Qi i‘^susy) 

= Diag(106.2,100.0, 81.4) TeV, 

(7.30) 

Tfl^c {rflsusy) 

= Diag(22.3, 22.3,49.0) TeV, 

(7.31) 

Tfl^C (yTflg^gy^ 

= Diag(138.0,132.1,120.1) TeV, 

(7.32) 

^Li i^susy^ 

= Diag(131.4,128.3, 84.4) TeV, 

(7.33) 


= Diag(172.0,164.4, 76.5) TeV, 

(7.34) 

ld>{jyisusy) 

= 91.6 TeV, 

(7.35) 

TflA (jT^susy) 

= 507.3 TeV, 

(7.36) 


Values of gauge and Yukawa couplings at Mqut: 

gf^ssMf^McuT) = (0.68367, 0.67083, 0.66722), 

Qgut = 0.71339, 

{KKK){Mgut) = (0.51487,0.03185,0.03132), 

(Ae, \^){Mgut) = (1.47 X 10-^ 1.86 x 10-^ 1.86 x lO’^), 
{KAdAe){MGUT) = (3 X 10-®, 9 X 10-^ 9 X 10-6), 

GUT state masses: 

rriT = 7.1 X 10^6 GeV, 

= 5.1xlO^^GeV, 
my = 8.8 X 10^6 GeV, 

Other parameters and observables: 

= 24.7, 

rp(p+ ^ K+v) = 3.9 X 10^^ yrs ( for ^ 

gmin ^ ^ 

corresponding to a vacuum lifetime of 

gmin 

Tvacuum = , ^ ,3 - ~ lO^^^ yrs . (7.48) 

+0 

'"susy ^universe 

For completeness, we also present the running of the A-terms in Fig. 5 and of the soft scalar 
masses in Fig. 6. 

Let us describe briefly the main features of this point. Most sfermions he between 50 
and 170 TeV, with however the ur and cr as light as 22.3 TeV. The gauginos are typically 


(7.42) 

(7.43) 

(7.44) 


— 02 — 03 — 0 ) , 


(7.45) 

(7.46) 

(7.47) 


(7.37) 

(7.38) 

(7.39) 

(7.40) 

(7.41) 
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Figure 5: Running A-terms and gaugino masses for the parameter space point described 
in the text (magenta dot in Fig. 4), obtained using the 1-loop MSSM RGBs. Top and 
bottom left hgures: 1®*, 2 "’^ and 3^*^ generation A-terms. The blue, yellow and green lines 
represent A„., and A^,., respectively. Bottom right hgure: the blue, yellow and green 
lines represent Mi, M 2 and M 3 , respectively. 

lighter than the sfermions and the higgsinos, and the lightest supersymmetric particle 
(LSP) is the bino, with mass Mi = 16.9 TeV (the ur and the cr are the co-NLSPs, i.e. the 
next-to-lightest supersymmetric particles). The superpartners are therefore out of reach of 
present and next-generation colliders, and supersymmetric contributions to flavour physics 
observables are strongly suppressed, but proton decay will be easily accessible at future 
large detectors. The Higgs sector is far into the decoupling regime, with a standard-like 
lightest Higgs boson and all non-standard Higgs bosons around tua ~ 500 TeV. Finally, 
there is no suitable dark matter particle in the observable sector, as the relic density of 
a 17 TeV bino by far exceeds VtcDuh'^ — 0.12 [47]. If R-parity is conserved, a natural 
candidate is a gravitino in the GeV range (or lower), to which the bino NLSP would 
decay without affecting Big Bang nucleosynthesis. The cold dark matter density would 
then be in the form of gravitinos coming from NLSP decays [48] (giving a contribution 
AGg = fl^niQ/Mi, from which the upper bound tuq < 1 GeV follows) and from thermal 
production during reheating [49, 50]. In order to avoid gravitino overproduction from the 
latter process, the reheating temperature should lie in the TeV range, or lower. In the case 
of R-parity violation (which may be invoked to generate neutrino masses, as an alternative 
to the seesaw mechanism), the only possible dark matter candidate within this model is 
again the gravitino. Since it decays into a photon and a neutrino, acting as a source of 
monochromatic photons, extragalactic gamma ray constraints put a bound on its mass. 
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generation soft masses squared 
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Figure 6: Running soft scalar masses for the parameter space point described in the text 
(magenta dot in Fig. 4), obtained using the 1-loop MSSM RGBs. Top and bottom left 
figures: 1^*, 2”“^ and 3'”'^ generation soft sfermion masses. The blue, yellow, green, red and 
purple lines represent m% , mL, m|c, and respectively. Bottom right hgure: soft 

Higgs boson masses. The blue line represents and the yellow line 

For bilinear R-parity violation, one obtains tuq < 1 GeV [51]. Here, as in the R-parity 
conserving case, the maximal reheating temperature is around the TeV scale. 

One may wonder whether other points exist in the parameter space of the minimal 
renormalizable SU(5) model with lighter superpartners than in the above example. While 
performing an extensive scan would probably reveal the existence of such points, it appears 
difficult to lower signihcantly the average superpartner mass scale below a few tens of TeV 
(some supersymmetric particles may however be accidentally lighter). The reason for this 
is the proton lifetime, which in the limit where all superpartners have the same mass rususy 
scales as Tp oc rn\rnl^gy oc rn\l/^ (where in the last step the approximate constraint (3.19) 
from gauge coupling unihcation was used). The actual proton decay constraint depends on 
the individual superpartner masses^^, and does not prevent some of them to be much smaller 
than the naive lower bound on rUsusy, but it seems to be difficult to reconcile a signihcantly 
lighter superpartner spectrum with proton decay and all other constraints. In fact gauge 
coupling unihcation is the main obstacle here, which makes it hard to accommodate a low 
superpartner mass scale with a large colour triplet mass [7]. 

What one may still try is to split the spectrum, making just a part of it light. This is 
not easy though, because of the threshold corrections to gauge couplings at Mqut and to 
the bottom and Higgs masses at rUsusy For example, making the higgsino light decreases 

addition, non-zero values of the high-energy phases in the Yukawa matrices tend to increase the 
proton lifetime. 
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the mass of the colour triplet, see Eq. (3.12), and therefore shortens the proton lifetime. 
The gluino and/or the cannot be much lighter than the (whose mass is of order rrisusy) 
either, since this would suppress the loop function in Eq. (3.26) and require an enormous 
value of Xb to reproduce the observed mb/rrir ratio. This large value of Xb would in turn 
imply a large negative threshold correction to the Higgs quartic coupling, see Eq. (3.38), 
making it impossible to £t the measured Higgs mass for moderate values of nisusy All 
this motivates the choice of parameter intervals made in our scan, namely 'mQ_^ ^{nisusy), 
mga {rrisusy) and M 3 {Mgut) in fhe range [0.5 rrisusy, 2 rrisusy]- The only remaining possibility 
is to have the hrst two sfermion generations lighter than the third one, an opposite situation 
to the one considered e.g. in Ref. [7]. The investigation of this case would require the 
use of 2-loop RGEs for soft terms, as a strong mass hierarchy between different sfermion 
generations enhances the effect of 2-loop running and may give rise to tachyons [52]. Such 
a study is beyond the scope of this paper, and we leave it for future work. 

8 Conclusions 

If one excepts the issue of neutrino masses, the minimal renormalizable supersymmetric 
SU(5) model suffers from two main problems. The hrst one concerns the predictions for 
charged fermion masses. Although the GUT-scale equality of down-type quark and charged 
lepton Yukawa couplings leads to a qualitatively successful prediction for the bottom to tan 
mass ratio, it quantitatively differs from the measured value by some 20-30%, while the 
discrepancy is much larger for the hrst and second generations. Second, proton decay is too 
fast by a factor 10^ or so for typical soft terms in the TeV range, essentially because gauge 
coupling unification requires a relatively light colour triplet. Each of the two problems 
has been addressed separately in the literature: large supersymmetric threshold corrections 
were used to modify the fermion mass predictions (see for example Refs. [11, 12]); specific 
havour structures of the soft terms [3] and a heavy superpartner spectrum [53] were invoked 
to increase the proton lifetime. 

In this paper, we performed a complete analysis of the minimal renormalizable super- 
symmetric SU(5) model, taking into account all relevant phenomenological and theoretical 
constraints: charged fermion masses, the proton lifetime, gauge coupling unihcation, the 
Higgs mass, experimental lower bounds on superpartner masses and havour constraints, 
metastability of the vacuum and perturbativity of the model. We showed that the model is 
still alive, and that the allowed region of the parameter space spreads over a large domain 
with rrisusy = ranging from around 50 TeV to 10^ TeV and tan (3 from approxi¬ 

mately 2 to 5. Viable points were also found outside this region of the parameter space, 
but they are less frequent. Particularly interesting are the ones featuring some (relatively) 
light superpartners. We studied in greater detail one such point in which the lightest 
supersymmetric particle is a bino with mass 16.9 TeV. 

A generic feature of the model is the metastability of the vacuum, which is not a concern 
since its lifetime is typically much larger than the age of the universe. A consequence of 
the heavy supertpartner spectrum is that the observable sector does not contain a suitable 
dark matter candidate; however, a gravitino with a mass around the GeV scale or below 
can play the role of cold dark matter, both in the presence and in the absence of R-parity. 

The present analysis is a good starting point for further research. One of the limitations 
is that we used the central values of the input SM parameters; taking into account the exper¬ 
imental and theoretical uncertainties will extend the viable region of the parameter space. 


30 



We presented in a compact way all the ingredients that are needed for a more extensive 
scan of the parameter space of the model. What we have shown is that the minimal renor- 
malizable supersymmetric SU(5) model is indeed a viable extension of the Standard Model, 
in the sense that its parameter space contains points that satisfy all relevant phenomeno¬ 
logical and theoretical constraints. The price to pay is that the superpartner spectrum is 
heavy, typically in the 100 TeV region or above, although we were able to hnd particular 
cases in which some of the supersymmetric particles can be as light as 17 TeV, and it is 
not excluded that viable points with even more split spectra exist in some corners of the 
parameter space. 

We did not address in this paper the issue of neutrino masses. One possibility is to add 
Standard Model singlets (right-handed neutrinos) to the model and to generate neutrino 
masses via the seesaw mechanism. Another possibility is to allow either bilinear or trilinear 
R-parity violating terms. Notice that integrating out higgsinos and colour triplets with 
R-parity violating couplings induces corrections to the Yukawa couplings [21], which opens 
the possibility that the SU(5) fermion mass relations are cured by the combined effect of 
R-parity violation and of supersymmetric threshold corrections. 

The possibilities to test experimentally the minimal renormalizable supersymmetric 
SU(5) model are limited. Barring possible isolated points of the parameter space with light 
remnants, superpartners are too heavy to be detected at the next generation of colliders 
or to give sizable contributions to flavour-violating observables. Proton decay could be 
around the corner, but the allowed parameter space is too vast to give a dehnite prediction 
for the proton lifetime. It is, however, possible to rule out the model (or at least most of 
its allowed parameter space), either by discovering a light superpartner spectrum or, if the 
Higgs sector is light enough to be accessible and studied at the next generation of colliders, 
by measuring a relatively large value of tan/?. 
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A Procedure for solving the soft term RGEs 

In this appendix, we derive the expressions for the superpartner spectrum of the minimal 
renormalizable supersymmetric SU(5) model that we have used in our analysis. These 
formulae account for the splitting of sfermion masses within the same SU(5) multiplet due 
to renormalization group running. They are (approximate) semi-analytic solutions to the 
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1-loop MSSM RGEs for the soft terms obtained under certain assumptions (namely, the 
hierarchy among Yukawa couplings allows us to neglect some terms in the RGEs and to 
solve them in a sequential manner, as explained below). We have checked that they are 
accurate up to the few percent level by running the full set of RGEs. The advantage of 
these approximate solutions is that the low-energy soft terms can be written as quadratic 
functions of the initial (GUT-scale) values of the soft parameters, which makes it possible 
to perform a scan over the parameter space of the model without having to solve the full 
set of RGEs for each point. Furthermore, it is possible to scan directly over the low-energy 
values of the soft terms (rather than the GUT-scale ones), since these formulae implicitly 
respect the SU(5) boundary conditions. A similar approach has already been employed in 
the past, e.g. in Ref. [26] for mSUGRA and in Ref. [27] for SU(5). Here we generalize this 
procedure to generation-dependent soft terms with SU(5) boundary conditions. 

A.l The procedure 

In order to perform a scan over the parameter space of the minimal renormalizable su¬ 
persymmetric SU(5) model, one has to solve a system of entangled differential equations 
describing the running of gauge, Yukawa and Higgs quartic couplings, as well as of the soft 
terms, which involves a large number of parameters (only in the supersymmetry breaking 
sector, there are upon SU(5) unihcation 15 parameters, namely 6 A-terms, 1 gaugino mass 
and 8 soft scalar masses). In addition, these parameters are subject to various phenomeno¬ 
logical constraints defined at different scales. This makes it hard to solve the problem by 
brute force, and motivates the use of approximate semi-analytical solutions to the soft term 
RGEs. 

Let us first try to circumvent the issue by solving the RGEs in steps. We first integrate 
numerically the system of RGEs for gauge, Yukawa and Higgs quartic couplings at the 2- 
loop level. Note that these couplings depend on the soft terms only through supersymmetric 
threshold corrections, so we can use this fact to determine the values of the relevant soft 
terms at the matching scale mg^sy Then we can solve the 1-loop RGEs for the soft terms 
- first gaugino masses, then A-terms and finally soft scalar masses. By working them 
out in this particular order, we are able to use in each step the knowledge of the running 
parameters that we have computed in the previous steps. Unfortunately this method is not 
efficient since one must solve the RGEs for the soft terms every time one changes the input 
(GUT-scale) values of the soft parameters. What we would like is to numerically solve the 
soft term RGEs for symbolic input values. The computational procedure described below 
does precisely that: when the input values of the soft parameters are changed, only the last 
step (out of 3) has to be re-run. 

Step 1 

1. Ghoose the matching scales nigusy and Mguti the ratio of the MSSM Higgs vevs tan [3 
and the masses of the heavy GUT states rriT and my within the perturbative regime 
rriT < rny -C Mpianck (note that rriT and my affect the running of the parameters 
only through the high-scale threshold corrections to the Yukawa couplings, so their 
influence is very mild). 

2. Run the gauge, charged lepton and up-type quark Yukawa couplings as well as the 
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Higgs quartic coupling from their measured values at the weak scale up to the match¬ 
ing scale rususy using the 2-loop SM RGBs. 

3. Apply the matching conditions (3.5)-(3.7), (3.21) and (3.23) to the gauge and Yukawa 
couplings at the scale msusy and use the 2-loop MSSM RGBs to run them up to the 
scale Mgut- 

4. Gompute the down-type quark Yukawa couplings at Mgut from the values of the 
charged lepton Yukawa couplings using the SU(5) mass relation (3.25) and the GUT 
threshold corrections (3.35). Run them down to the scale msusy using the 2-loop 
MSSM RGBs (the supersymmetric threshold corrections to the leptonic Yukawa cou¬ 
plings are neglected in this procedure). 

Step 2 

5. Solve the 1-loop MSSM RGBs for gaugino masses, taking into account the unihcation 
condition (3.46). 

6 . Assume the hierarchy^^ A* > A?, > A,- > Ac > As > > Ag and the 

unihcation of leptonic and down-type quark A-terms at the GUT scale to solve the 
RGBs for A-terms in a sequential way (see details in Subsection A.3. 2). As a result, 
express the running A-terms as linear functions of the gangino masses and of their 
valnes at the scale mgusy (which are more convenient inpnt parameters than their 
valnes at the GUT scale, since phenomenological constraints on A-terms apply at 
the scale msusy), with coefficients depending on the running qnantities determined in 
Step 1. In this way one does not need to rnn the A-terms every time one changes 
their initial values. 

7. Solve the RGBs for the soft scalar masses in a similar manner. 

Step 3 

8 . Ghoose the values of the input soft parameters. 

9. Make sure that no sfermion becomes tachyonic. Gheck that the values of the A- 
terms and soft scalar masses at msusy reprodnce the observed values of the down-type 
quark masses and of the Higgs mass. Verify that all experimental lower bonnds on 
superpartner masses and flavonr constraints are satisfied. 

10. Gheck vacunm (meta)stability constraints and the proton lifetime. 

^^This assumption is based on the observed hierarchy of fermion masses, taking into account the SU(5) 
boundary condition Ajj = and the supersymmetric threshold corrections to down-type quark masses, 
and is better justified in the low to moderate tand regime. For instance, the relation Xs{Mgut) = 
\^{Mgut) implies Xs{m) > X^{m) for msusy < m < Mgut, while in the absence of this relation and 
of supersymmetric threshold corrections to mg one would have Xs(jn) < A^(m) over most of the range 
msusy <m < Mgut- 
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A.2 Input values 


All input values of the SM parameters needed for the running were taken either from 
Ref. [54] or from Ref. [55]: 


^Planck 

V 

9i{mt) 

{mc,ms,m^){mz) 

{mu,md,me){mz) 

A(mt) 


2.4 X lO^^GeV, 

174.10362 GeV, 

(91.1876, 173.10, 125.66) GeV, 
(0.46167, 0.64822, 1.1666), 
(0.93558, 0.0156, 0.0100), 

(619, 55, 103) MeV, 

(1.27, 2.9, 0.5) MeV, 

0.25420, 


(A.l) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 


where S'!, 2 , 3 (^ 4 ), and \{mt) are evaluated at NNLO (the NNNLO pure QGD 

contribution is also included in htirrit) [54]), and hiijnz) = mi{mz)/v for the hrst and 
second generations of fermions. 


A.3 Approximate expressions for the soft terms 

In order to be able to approximately solve the 1-loop RGBs for soft terms, we shall neglect 
all mixings. This means that in addition to assuming that the sfermion soft terms are 
aligned with fermion masses at the GUT scale (see Subsection 2.1), we shall neglect the 
effects of the GKM matrix in the RGBs. This may not be fully justihed, as the impact of 
VcKM on the running of the hrst two generation parameters can be signihcant. However, 
this effect is suppressed either by small Yukawa couplings or by small GKM angles and is 
therefore never numerically important^'^, so we shall set Vckm = 1 in the RGBs and omit 
the subleading Yukawa contributions. 

With these assumptions, we can derive approximate semi-analytic solutions to the 1- 
loop RGBs for the soft terms. 


A.3.1 Gaugino masses 

The 1-loop RGBs for gaugino masses: 


d ,, , da^ Mi 

—Mi = ^- aiMi ~ ^ , 

dt dt (Xi 


(A.9) 


(where the last equality is only approximate because we run gauge couplings at 2 loops) 
has the following simple solution respecting SU(5) symmetry: 


Mi[m) = Mi/2e Jinm ^ ' 


- ^ 1/2 y-yyy-^ - Miirrig) ^ ^ ~ rrig ^ Z~TT7 -1 • (A.IO) 


cxiiMcup) 


ai[mg) 


Oisirrig) ai^McuT) 


This solution becomes exact when the 1-loop RGBs for gauge couplings are used. 


^^Obviously, this statement does not apply to RG-induced flavour-violating soft terms. These are not a 
concern, however, since the superpartner spectrum is heavy and the first two generation squarks are almost 
degenerate in mass. 
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A.3.2 A-terms 


For small to moderate tan/?, one can take advantage of the hierarchy among Yukawa 
couplings (see Subsection A.l, Step 2) to simplify the 1-loop RGBs for A-terms and to 
solve them in a sequential way^®. Namely, one can write the A-term RGBs as: 


d , _ An 

Jt " ~ (4^ 




L P 


i=l 


+ 


^r. 


(47r)" 


p^n 2=1 


, (A.ll) 


where the indices n and p run over the ordered values {t, b, r, c, s, p, u, d, e} (such that e.g. 
p < c means p = t,b or t) and the coefficients J”, are collected in Table 1. 


p 

t 

b 

r 

c 

s 

/i 

u 

d 

e 

i 

3 

2 

1 

p 

v 

18 

1* 


3* 



3* 



Kl 

16 

3 

3 

13 

15 

Jb 

1 

18 

1* 


3* 

1* 


3* 

1* 

Kt 

3 

3 

7 

15 



3 

12 


3* 

1* 


3* 

1* 

KJ 


3 

9 

5 

^p 

3 



18 

1* 


3* 



Kf 

16 

3 

3 

13 

15 



3 

1 

1 

18 

1* 


3* 

1* 

Kf 

16 

3 

3 

7 

15 

jp 

■jp 


3 

1 


3 

12 


3* 

1* 

AT 


3 

h 
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JU 

3 



3 



18 

1* 


Kf 

16 

3. 

3 

13 

15 

jd- 

•Jp 


3 

1 


3 

1 

1 

18 

1* 

Kf 

16 

3 

3 

7 

15 

je 


3 

1 


3 

1 


3 

12 

Kf 


3 

9 

5 


Table 1: The coefficients appearing in the 1-loop RGBs for A-terms (A.ll). 

A superscript * on a coefficient in Table 1 indicates that we neglect the correspond¬ 
ing term 2AnApAp/(47r)^ in the RHS of Bq. (A.ll), consistently with the hierarchy of 
Yukawa couplings^®. The evolution of each A-term is then (approximately) described by a 
differential equation of the form 

4 am = /3„yi)A„(i)+7^W. (A.12) 

where the coefficients /3:^(t) and 7^(t) depend on the gauge and Yukawa couplings, and 
in addition depends linearly on the gaugino masses and on the A-terms Am(t) with 

m < n. This makes it possible to solve the RGBs (A.12) sequentially, from n = t to n = e. 
The solutions can be written as: 


An{t) = 


+ 




dt' 


no 




(A.13) 


where, by using the already obtained expressions for the A^’s with m < n, the integrals 
of the coefficients (d^{t) and 7,^(t) can be computed numerically after having solved the 

case the Yukawa-dependent terms in the RGEs do not follow the hierarchy assumed for the Yukawa 
couplings themselves, one can if necessary improve the accuracy of the solutions by iterating the procedure 
described below. 

^®In addition, one does not need to include the terms that are suppressed by small Yukawa couplings in 
the RGEs (A.ll), as their effect is smaller than the precision of the 1-loop approximation. We nevertheless 
give the corresponding coefficients in Table 1 for completeness, allowing for the possibility of unusually 
large first or second generation A-terms that would make some of these terms relevant. 
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MSSM RGBs for the gauge and Yukawa couplings. As a result, one obtains the running 
A-terms at an arbitrary scale as linear combinations of the 7 SU(5) soft parameters M 1 / 2 , 
Al^, Af {i = 1,2,3), with numerical coefficients depending on the choice of rngusy and 
tan/? (and very mildly on Mgut, 'tut and my). In practice, it will prove convenient for 
the exploration of the parameter space of the model to trade these GUT-scale parameters 
for low-energy ones, so as to express the An{t) as a function of M^{msusy)i Ad^i^rUsusy) and 
AuX^susy) (* = h 2, 3). Note that if the SU(5) boundary conditions at Mgut had not been 
imposed, the An{t) would depend on 12 initial parameters, namely Mi 2 , 3 , A„., Ad- and Af,-. 


A.3.3 Soft scalar masses 


One can apply a similar procedure to the 1-loop RGBs for soft scalar masses. Let us hrst 
dehne the following variables: 


> 

(A.14) 

’ 

(A.15) 

2 1 2 1 2 

mi,+ mgc + rriH^ . 

(A.16) 


These combinations of masses appear on the RHS of the RGBs for soft scalar masses, and 
obey 1-loop RGBs of the form 


1 

dt 


jn\2 


(dvr) 


+ 


(dvr) 


p^n P 




i=l 


(A.17) 


where again the indices n and p run over the ordered values {t, h, r, c, s, p, u, d, e} and the 
coefficients L” are collected in Table 2, while =‘^Jp and A” = dAf. 


n 

t 

b 

r 

c 

s 

h 

u 

d 

e 

Tn 

12 

12 

8 

12 

12 

8 

12 

12 

8 


Table 2 : The coefficients L” appearing in the 1-loop RGBs for (A.17). 

Neglecting^^ the terms LpApEp/(d 7 r)^ (p 7 ^ n) on the RHS of Bq. (A.17) when the coef- 
hcients J” are marked with a superscript * in Table 1, one can solve the RGBs for the 
E„’s sequentially from n = t to n = e, as we did for the A-term RGBs. Indeed, in this 
approximation Bq. (A.17) can be written in the form: 




(A.18) 


where is proportional to and 7 ^(t) is a linear combination of gf(t)Mf(t), A^it) 

and A^(f)Em(t) with m < n. The solution reads: 

E„(i) = E„(i„) eh J‘dt' , (A.19) 


where the integrals involve only already known quantities, since the RGBs for the E^’s are 
solved in order of increasing m. Using the previously obtained expressions for the running A- 
terms and gaugino masses, the combinations of masses E„(t) are then expressed as quadratic 

^^The terms suppressed by Yukawa couplings of the first and second generations can also be dropped, 
as their effect is smaller than the precision of the 1-loop approximation. 
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functions of the GUT-scale soft scalar masses mio^, mg., mnS^GUT)) '^H^iMcuT) and of 
the nisusy values of M 3 , Au,c,t and Ad,s,b, with t-dependent coefficients that can be computed 
numerically using the known solutions to the 2-loop MSSM RGBs for gauge and Yukawa 
couplings. 

Another combination of masses that appears in the RGBs for soft scalar masses is 


^ = 


— 


m 


Hu 


m 


Hd 


+ Tr 


m 


Q 


m, 


2m^ 


+ ml 


+ m| 


= ^Hu -^Hd + T. hi. “ ~ , 

i=l 

for which the 1-loop RGB takes the simple form: 

d 66 «! ]jMSSM g 

—5 = — — 5 = -ctih ~ —-. 

dt 5 dvr 27r dt ai 


(A. 20 ) 


(A. 21 ) 


The last equality is only approximate because we run gauge couplings at 2 loops. Working 
in the 1-loop approximation, one can solve Bq. (A.21) straightforwardly: 


pin A 

5'(m) = S{^MqU'J'^ € Jlnm 

ai(m) 


— S{Mgut) 


a gut 


dt Qi (t) /{27 v) 


I^Hui^GUr) —Hfi^jiMGUT) 


- aiijn) 
J a gut 


(A. 22 ) 


where we have used the fact that S{Mgut) = (M^j/t) — hi%^{Mgut) due to the SU(5) 

boundary conditions on soft scalar masses. 

We now have all the ingredients needed to write the solutions to the 1-loop RGBs for 
soft scalar masses as: 


ml{m) = mliMGUT) + ^P“/s,,(m) + ^P“/A„(m) - ^Q"/M,(m) + R^Is{m), (A.23) 

n n i=l 

where a = Hu, Hd,Qi,u^,e1,Li,d1, the numerical coefficients P", Qf and P“ are given in 
Table 3, 


a 

pa 

pa 

pa 

r,Ai,e 

Qs 


Q? 

P" 

Hu 

6 




6 

6 

5 

+i 

Hd 


6 

2 


6 

6 

5 

3 

5 


a 

pa 

^Ui 

pa 

^d. 

pa 

Q 3 

Qf 

Qf 

P" 

Q^ 

2 

2 


32 

3 

6 

2 

15 


u1 

4 



I 2 

3 


32 

15 

4 

e? 



4 



24 

5 

+ 5 

R 



2 


6 

6 

5 

3 

5 

dt 


4 


32 

3 


8 

15 

+i 


Table 3: The non-vanishing coefficients P", Qf and P" appearing in Bq. (A.23). 
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and the integrals /s„, Ia„: I Mi and Is are defined by: 


lAn{m) 

lMi{.m) 

Is{m) 


(dvr) 2 
1 

(dvr)^ 

1 


^In m 

hnmcuT 

^In m 

hnmauT 

/•In m 


dtAl{t ), 


Ml,, 


(dvr)^ 

1 

^dyr^^ JinTtiQU']^ 


In mauT 

/•In m 


dtg,{t)M^{t) ~ 

S{Mgut) 


dtg^{t)S{t) ~ 


21)Mssm 


^ Oil {Mgut) 
ai{m) 
Oii{MGUT) 



(A.2d) 

(A.25) 

(A.26) 

(A.27) 


Plngging the semi-analytic expressions for S„(t) and Anit) into Eqs. (A.2d) and (A.25), 
one can express the rnnning soft scalar masses as qnadratic fnnctions of the GUT-scale soft 
parameters mio,, mg., mH^{MGUT), rnH^{MGUT) and of the rrisusy valnes of Mg, Au^c,t and 
Ad,s,b, with m-dependent coefficients that can be compnted nnmerically. 

One can further simplify Eq. (A.23) by writing /s„ as a linear combination of the 
differences Ep(m) — T,p[Mgut) and of the integrals Ia^ and I Mi- Indeed, satisfies a 
1-loop RGE of the form: 

dE„(() = ;3j(«)E„(()+7„y«), (A.28) 


with /3n{t) equal to 12A^(t)/(d7r)^ for squarks and to 8A^(f)/(d7r)^ for sleptons. From this 
and from the definition (A.2d) of /s„, one immediately deduces that 


^Se, {m) 


~ ^Ui,di{MGUT) 

12 

Eei(m) — 'Eei{MGUT) _ 1 
8 8 


plnm 

^In m 
/In Mqujp 


(A.29) 

(A.30) 


where the integrals on the RHS can be expressed in terms of Js^, Iap and I Mi, so that 
eventually all /s„’s can be written as linear combinations of T,p{m) — T,p{Mgut), Iap and 
I Mi- This leads to the following expressions for the running soft scalar masses: 


ml{m) = ml{MGUT) + T.T:{^n{m)-J:n{MGUT)) 

n 

+ E KTAm) + E V^MSm) + RVsim ), (A.31) 

n 1=1 


with the numerical coefficients T", G", (resp. i?“) given in Table d (resp. Table 3). 

The quantities S„(m), /A„(m), /^^.(m) and /^(m) on the RHS of Eq. (A.31) - hence the 
running soft scalar masses m^(m) - are quadratic functions of chosen initial parameters, 
with coefficients depending on the scale m, tan (3 and nisusy (and very weakly on Mgut, 
and my)- In the above derivation, the initial parameters were taken to be the GUT-scale 
soft masses mio,, mg., uih^^Mgut), mH^McuT) and the rrisusy values of Mg, and 

Ad,s,b- Alternatively, one can trade the 8 GUT-scale masses mio,, mg., rriH^^MGUT) and 
rriH^iMGUT) for 8 low-energy soft scalar masses by inverting 8 of the 17 equations (A.31), so 
as to express all running soft parameters (soft scalar masses, gaugino masses and A-terms) 
as functions of the rrisusy value of 15 of them. 
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23 
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3 
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3 

32 
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54 

59 

96 
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Table 4: The non-vanishing coefficients T", and 14“ appearing in Eq. (A.31). 


B Proton lifetime computation 


In this appendix, we outline the computation of the higgsino-mediated (D=5) proton decay 
rate [56, 57, 33, 58] in the minimal renormalizable supersymmetric SU(5) model. Integrating 
out the heavy colour triplet components of the bn and bu superhelds generates the following 
D=5 operators: 


1 

2 


^ijkl 

TUt 


QiQ jQkLl 



1 

2 rriT 




h.c., 

6»2 


(B.l) 


where the parameters ntjki and can be expressed in terms of the SU{5) Yukawa cou¬ 
plings and A|- dehned in Eq. (2.2): 


K^ijki — Ab’Affc, — AjJ^A^j, (B.2) 

and the contraction of gauge indices is understood. Since colour invariance implies Kmi = 0 
and = 0, the dominant proton decay modes arising from the above operators involve 
a kaon, and in practice p —)■ dominates. The corresponding amplitude is obtained by 

“dressing” the D=5 operators of Eq. (B.l) with gaugino/higgsino loops. Over the region of 
the parameter space considered in this paper, the dominant contribution comes from the 
wino dressing of the QQQL operator^®, and the decay rate reads: 
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- 

— I'^kj 


A10a5 


T{p 


ly = 


mp ( m^+ 


327r 


mi 


E 

I 




1121 


+ C 


12U 


(B.3) 


^®The gluino dressing of the QQQL operator can be neglected as the mass difference between the first 
two generations of squarks is small, and the charged higgsino dressing of the u'^u'^<Fe‘^ operator is significant 
only for large values of tan (3. Although we consider large A-terms, the left-right sfermion mixing remains 
small <C 1), hence the wino dressing of the operator is not relevant either. 
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where Cijki{fi) is the Wilson coefficient of the 4-ferniion operator Uidjdki'i and the hadronic 
parameter is defined by {K^\{uLidLj)dLk\p) — Wq^^PlUp, where Up is the proton 

spinor. The lattice computation of Ref. [59] gives = 2 GeV) = (0.111 ± 0.027) GeV^ 

and = 2 GeV) = (0.036 ± 0.014) GeV^, with statistical and systematic errors added 

in quadrature. In order to minimize the dependence of the proton decay rate on the 
renormalization scale p, the Cijkis are evaluated at the scale /i = 2 GeV rather than at 
fi = rrip. 

The Wilson coefficients Cijki are computed at the matching scale rUsusy at which su¬ 
persymmetric partners are integrated out by “dressing” the QQQL operators with loops 
containing a wino: 


a 


ijkl 




susy ) 




m,n,p 


^pmnl 

rriT 


Y, V:Ypy nw, dLg, ULm) + /(IV, Uln, Hi) , (B.4) 


in which V is the GKM matrix, the parameters Kpmni = f^pmni — i^nmpi are expressed in 
the super-GKM basis in which the up quarks are mass eigenstates, and the possible mis¬ 
alignment between the fermion and the sfermion mass eigenstate bases is neglected. In 
Eq. (B.4), the GKM matrix entries and the KpmnH are renormalized at the scale rngusy 
The loop function /(IV, a, b) is given by: 

/(lV,a,6) = ^m*/3(m^,mb,m|) , (B.5) 


where the function I^{x,y,z) has been dehned in Eq. (3.30). For equal sfermion masses 
(rua = mt = rrisusy), /(IV, a, 6) reads: 


/(IV, a, 6) 


02 1 
1- b 

AlTI ‘kklgusy 


( ^lusy 


9{.x) 


a/x (x — 1 — In x) 
(x — 1)2 


(B.6) 


The parameters hijki depend on the details of the Yukawa sector of the Grand Unihed 
Theory. In the minimal supersymmetric SU{5) model, they are given by (neglecting the 
running between Mqut and the triplet mass scale): 




A, 


.e 






(B.7) 


where the Yukawa couplings and the GKM matrix are evaluated at the scale m^, and the 
are high-energy phases satisfying the constraint 0i -l- 02 + 03 = 0. These parameters 
must then be evolved down to the scale rrisusy by solving the appropriate RGEs. Neglecting 
the Yukawa couplings, the running simply amounts to an overall rescaling: 


^ijklij^susy) -^SD ^ijklip^T) ) (B-S) 

where in the minimal renormalizable supersymmetric SU{5) model: 

^ _ / Q:i(m^At^^) y^ / a2{msusy) y^ f «2(»^3) VV asirrisusy) 

\ aiirriT) J V ) \a2{mT)) V «3("iT) 

Finally, the Wilson coefficients Cii 2 i{rnausy) and Ci 2 ii{.''n^susy) must be renormalized down 
to the scale p = 2 GeV before being inserted into Eq. (B.3): 

Cii 2 i{p = 2 GeV) = Cii 2 i{msusy ), 

Ci 2 ii{p = 2 GeV) = Ci 2 ii{msusy ), 
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(B.9) 
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where and A^^^ are renormalization factors given by (nsing the formnlae of Ref. [60]): 
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The fact that A]^'^ ^ A^^^ is dne to the RG-indnced mixing between 4-fermion operators 

with different flavonr indices; when this (nnmerically small) mixing is neglected, A]^^ = 
4121 — A 

^LD - ^LD- 


C UFB constraints 


In this appendix, we derive the constraints on the model parameters associated with the 
absence of directions along which the scalar potential is nnbonnded from below. In doing 
this we will follow Ref. [61]. 


C.l UFB-1 

The absence of this direction (where Hu = Hd to cancel the D-terms and all other vevs are 
vanishing) reqnires 


+ > 25, 

(C.l) 

which can be rewritten (for finite mz) as 


22 2 ~ 1 

mri — uitt > rriy --, 

- ^tan2/3 + l’ 

(C.2) 

where 


, ,2 (3 m| 

tan2/3-l 2 ’ 

(C.3) 

^ - ^Hu ) tan /3 m| tan (3 

taiA (3 — 1 tsiiA (3 + 1 

(C.4) 

Eq. (C.2) is eqnivalent to 


5 > 0, 

(C.5) 

which is antomatically satished by assnmption. As a cross-check. 

one can notice that ne- 

glecting mz (since we typically consider sitnations where mfj 3> ni|), the constraint (C.2) 

rednces to: 


2 \ 2 
nriTT > ttItt , 
tid — riu^ 

(C.6) 


which is again antomatic since |/ip > 0. 
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C.2 UFB-2 

We now allow nonzero vevs for and i> only: 


r = + M^)h] + (mi. + 

(«l - Hi - 


9 “ + sl 


In the limit rriz —)■ 0, one can recast Eq. (C.7) into a sum of three positive terms: 
- m]j^ (ETrftan/? - H^f 


V = 


tan^ /3 — 1 


, 2~2 , + 

+ miiy + —-- 


(C.7) 


{hi - Hi - 7=) , (C. 8 ) 


hence there cannot be any minimum lower than the electroweak vacuum E = 0. For finite 
ruz only the first term in Eq. (C. 8 ) is modified, but this does not affect the minimization 
of V with respect to h: 


z/ 


u 



4m| \ 
9'^+ 92 ) 


0 . 


(C.9) 


The solution 0 = 0 brings us back to the UFB-1 case, so we are only interested in the non¬ 
trivial h 7 ^ 0 solution, which is however consistent only for > ifj. One can check both 
analytically (using for example the command Reduce in Mathematica) and numerically 
that in this case there is an UFB direction only if 


171% — m% < m% 


(C.IO) 


(plus some extra constraints), or 

(C.ll) 
(C.12) 

(plus some extra constraints). Since none of these conditions is satisfied in our analysis of 
the minimal renormalizable supersymmetric SU(5) model, we can conclude that there are 
no UFB-2 directions. 

Naively one could think that this conclusion may change if instead of h 7 ^ 0 one would 
allow a nonzero e = e'^, because it contributes with the opposite sign to the D-term. This 
is however not the case: while the non-trivial solution e = 7 ^ 0 is now consistent for 

Hi < Hj, there is also a sign change in Eq. (C.9). Hence also in this case there are no 
dangerous UFB-2 direction. 

C.3 UFB-3 

One can avoid the UFB-3 direction (with nonzero vevs for Hu, v and d = d^) if the 
conditions 

(C-13) 

are satisfied for any of the three slepton generations. 


TTIh^ — fTT' 


2 \ 2 

> m 


Hu - 


Z ) 

2 


2 rriz tan^ /S — 1 
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D Decays into CCB vacua 

In this appendix, we describe a method to estimate the lifetime of the electroweak vacuum 
in the presence of charge and colour breaking (CCB) vacua induced by large A-terms from 
any sector and generation. To this end, we will approximate the minimal value of the 
bounce action in Euclidean space [62, 63, 64]: 

fOQ 

S' = 27rM drr^ {Ckin + V) (0*(r)) (D.l) 

Jo 

for the solution of the equations of motion 

^"(r) + ^^'(r) = A , (D.2) 

r ocpi 

subject to the boundary conditions (describing helds sitting in the false vacuum at r = cx) 
and tending to a lower minimum at r = 0): 

0,(oo) = C, 0:(O) = O. (D.3) 

The general problem can be solved only numerically [16, 65, 17, 66, 67, 68, 69]. We will sim¬ 
plify it by minimizing the bounce action along a constant direction in held space [71, 72, 70], 
which effectively reduces the problem to a single differential equation. This approxima¬ 
tion provides an upper bound on the minimal value of the correct multi-held bounce ac¬ 
tion [73, 72]. This means that the points of the parameter space ruled out in this way are 
indeed characterized by a too short vacuum lifetime, but the method retains points that 
after a more careful and complete analysis would be rejected. 

What we have to do is therefore to minimize an Euclidean action of the form: 

S = dr “ ^30^ + ^0^) , (D.4) 

where the held 0(r) parametrizes a direction in held space, and check whether its minimal 
value is large enough for the lifetime of the electroweak vacuum to be larger than the age 
of the universe, namely 

^ > 400. (D.5) 

This is done by solving the equation of motion (D.2) and plugging the solution back into 
Eq. (D.4). A semi-analytical approximate solution for the action can be found in Ref. [74] 
(see also Ref. [75]): 

S=^S{k), (D.6) 

where 

R(0) ^ 45.4, (D.7) 

and 

. = Ag. (D.8) 

K < 0 indicates that the potential is unbounded from below, while k > 1/4 ensures that the 
minimum at 0 = 0 is absolutely stable. Hence we only have to check that the metastability 
constraint (D.5) is satished for 0 < k < 1/4. In this range of k values, one has: 

S{k) > R(0) ^ 45.4. (D.9) 
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D.l An example of CCB constraint 


As a concrete example, let us first consider the vacuum decay to CCB minima induced by a 
large A-term in the down squark sector [61]; then we will generalize the discussion to other 
A-terms. The relevant Lagrangian is: 


kin 

V 


= K{r)\^ + \K{rW + \d'{r)\^ + IWWP + |e'(r)p + |(e7(r)p , 
= 

+ (^Aditdhd - XdlJ^*d‘'dhl - Bhuhd + h.c.) 

+ (l^nP + 

+ ^ + \K\‘ + l|dT + IW? - |8p + 2|8f 

+ f i\h,? - - MT -18|A + f (MT - M'lA. 


(D.IO) 


where hu,d, d (d^) and e (e'^) stand for the neutral components of the MSSM Higgs doublets, 
the scalar partners of the left(right)-handed down-type quarks and of the left(right)-handed 
charged leptons, respectively. Generation indices are omitted, and all mixing angles are 
neglected. Finally, the parameters p and B are determined by electroweak symmetry 
breaking, i.e. in the limit mz —t 0 : 


Ihl 


2 


B 


tan^ (3 


J^u 


tan^ fd — 1 


(m 


Hd 


— m 


H, 


) tan B 


tan^ jB 


(D.ll) 

(D,12) 


From now on, all quantities are assumed to be real. This is equivalent to assuming that the 
relative phase between and p equals (modulo tt) the relative phase between Ad and B 
(a condition which may be realized in some scenarios of supersymmetry breaking), because 
in this special case all phases can be absorbed in the helds involved. 

For all CCB constraints except the one associated with At (which will be discussed 
separately) and the one associated with Ah in the large tan/3 regime (which will not be 
considered because, in the minimal renormalizable supersymmetric SU(5) model, large 
values of tan jB are excluded by a combination of constraints, see Section 7) one has A^ ^ g <C 
5 f ^2 3 ) therefore one can impose the constraint of vanishing SU(3), SU(2) and U(l) D- 
terms [61]. 

As explained at the beginning of the appendix, we will search for solutions to the held 
equations of motion of the form: 


{hu, hd, d, cF, e, e'") (r) = (i7„, Hd, d, d, e, e) 0 (r), 


(D,13) 


with Hu, Hd, d, e constant, where we have in addition imposed the SU(3) D-term constraint 
(which also implies that d and point in the same SU(3) direction) 


cF = d 


(D.14) 


(the alternative sign choice d^ = —d would just hip the sign of M 3 , see Eq. (D.18) below, 
without ahecting the action (D. 6 ), which depends on M|), as well as the compatibility of 
the U(l) and SU(2) D-term constraints: 


(D,15) 
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(D.16) 


Finally, the vanishing of the SU(2) D-term: 

+ = 0 , 

will be considered as a constraint in the following. 

With these assumptions, the coefficients of the bounce action (D.4) can be expressed 


as: 

M| = {rriQ + mj) dP + {m\ + m-) 

+ + l^t) Hi + {ml^ + Ifip) Hi - , (D.17) 

M 3 = —2(f{AdHd — Xd^J-Hu) , (D.18) 

z = 2(2d^ + 2e^ + Hl + Hfj, (D.19) 

A = \ld^{£ + 2Hj). (D.20) 

Notice that we chose B > 0 (see Eq. (D.12); this amounts to £x the sign of HuHa). Without 
loss of generality, we can also choose 

Ad>0, >0, (D.21) 


while the signs oi jj, and Hd can in principle be arbitrary. For a given point in parameter 
space (for which the signs of Ad and /i are hxed), one needs to consider both signs of Hd- 

D.2 General CCB constraints 

Let us now generalize the previous discussion to all CCB constraints (except for the one 
associated with At). The coefficients of the bounce action (D.4) can be written compactly 


as: 

+ (ml + |/i|2) + (m\ + |/i|2) - 2BHH , (D.22) 

M3(^, M, id) = -2ij\A^H - , (D.23) 

zi'P, H, H, c,) = 2 (2V>' + + ^2 ^ ^ P_24) 

X{ij,H) = Xlijlij^ + 2Hl, (D.25) 

and the D-term constraint is given by: 

0 = -^/;2 + (3 - 2c)^2 + H^-H\ (D.26) 


where the 6 different possibilities are summarized in Table 5 (the example considered pre¬ 
viously was case 2). 





H 

H 


M? 

b 

ml 

ml 

Cz 

c 

1 

d (= d") 

V 

Hd 

Hu 

m| + 

ml 


ml 

Htu 

1 

1 

2 

d (= d") 

e (= e'") 

Hd 

Hu 

m| + 

ml + m| 


^k 

2 

2 

3 

u (= u'^) 

V 

Hu 

Hd 

m| 

ml 

ml 

Hlu 

^k 

1 

2 

4 

u (= u'^) 

e (= e'") 

Hu 

Hd 

+ m| 

ml + m? 

^k 

^k 

2 

1 

5 

e (= e*") 

u (= u'^) 

Hd 

Hu 

m| -|- m? 


^k 

ml 

2 

1 

6 

e (= e'") 

d (= d'^) 

Hd 

Hu 

m| -|- m| 

+ ml 

^k 

ml 

2 

2 
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Table 5: The six different cases of CCB constraints considered. The D-term constraints are 
given in parenthesis. Whenever d and d'^ {u and are involved, they point in the same 
811 ( 3)0 direction. 

For each case, there are several subcases to be checked according to which generation 
indices are carried by the V’ and V’ fields. For a given generation xjji (the 2 fields belong 
to the same generation), one must consider all generation choices for the il) fields that do 
not give a substantial contribution to the F-terms due to the size of the corresponding 
Yukawa couplings (otherwise there would be no CCB minimum lower than the electroweak 
vacuum). Namely, 

• when a single '0 field is involved (Cases 1 and 3), all generation indices should be 
considered; 

• when two 'ijj fields are involved (Cases 2, 4, 5 and 6), all combinations of two different 
generation indices should be considered; 

• the cases of two '0 fields with the same generation index should also be considered 
when the Yukawa coupling of the corresponding lepton is much smaller than Ap (i.e. 
one needs to consider only generations of 'ip^s whose fermionic partners are lighter 
than the fermionic partner of 'ip). 

In practice, it is enough to check generations of the 'ijj fields leading to the smallest value 
of M|, since the bounce action increases with M|. There are thus 16 possible cases: 
3 generations for case 1 with = {d,s,b), 3 generations for case 2 with = {d,s,b), 2 
generations for case 3 with = (u, c), 2 generations for case 4 with -ip = (u, c), 3 generations 
for case 5 with tp = (e, /i, r) and 3 generations for case 6 with xp = (e, p, r). 

Finally, starting from the semi-analytic approximation (D.6), it is easy to show that: 


^ = 4^(k) 




+ 2^2 ((3 + c, - 2c)tp‘^ + H^) + ((3 + c, - 2c)^^ + 2H‘^y 




> 4^(/s:)-y , 


(D.27) 


which is very large unless k is very small. Thus only very small values of k may lead 
to a short vacuum lifetime, and we can approximate S{k) ~ S'(O) ~ 45.4. To a good 
approximation, the quantity that we have to minimize is therefore 


5 = 




5(0) 


~w 


{ip,7p,H,H,c, c,) , 


(D.28) 


subject to the D-term constraint (D.26), and we must check whether 


qmin > ^ Q 

~ 45.4 ^ ’ 


(D.29) 


hereafter referred to as the metastability constraint. This condition must be satisfied for 
all six cases in Table 5 to be able to conclude (within our approximate method) that the 
lifetime of the electroweak vacuum is longer than the age of the universe. 

Notice that one needs to minimize S with respect to two variables only. Indeed, one of 
the variables ip,tp,H,H is fixed by the D-term constraint (D.26), and the quantity (D.28) 
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only depends on ratios of fields. The minimization is performed numerically. A point of 
the parameter space is admitted if the bounce action satishes Eq. (D.29) in all cases listed 
in Table 5. 

The D-term constraint (D.26) is most easily taken into account in the minimization of 
S by choosing a suitable parametrization of its general solution, for instance: 

D.2.1 Cases 1, 4, 5 (c = 1) 

For c = 1, the most general solution of Eq. (D.26) can be parametrized as: 


H = shH cosh.6 cos a , 
ip = H cosh 6 sin a , 

Ip = H sinh 9 , 


(D.30) 

(D.31) 

(D.32) 


with sh = ±1, 9 G [0, cx)[ and a G [ 0 , 7 r/ 2 ]. The actual value of 7 ^ 0 is irrelevant, as 
(D.28) only depends on ratios of helds. 

D.2.2 Cases 2, 3, 6 (c = 2) 

For c = 2, the most general solution of Eq. (D.26) can be parametrized as: 


H = shH cosh 9 , 

Ip = H sinh 9 sin a , 
'ip = H sinh 9 cos a , 


(D.33) 

(D.34) 

(D.35) 


with again sh = ±1, 9 E [0, cx)[ and a G [0, vr/2]. As in the previous case, the actual value 
of hf 7 ^ 0 is irrelevant. 


D.3 The special case of the stop 

In this case too we will make suitable approximations for estimating the minimal value 
of the bounce action (for more details on the subject, see Refs. [17, 76, 77, 78, 79]). As 
explained in Ref. [61], one cannot assume exact cancellation of the D-terms here, but we 
will nevertheless stick to the relatively good approximation of vanishing SU(3) D-terms, 
namely i = Keeping the contributions of all other helds, one can write: 

= \K{r)\^ + \K{r)\^ + 2\i\r)\\ 

+ mf) 

+ ~ — Bhuhd + h.c.^ 

+ (1^1^ + 




(D.36) 


Assuming again all quantities to be real, we take the same ansatz as for the other CCB 
constraints: 

{K,hd,i){r) = {Hu,Hd,t) (p{r), (D.37) 
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with Hu, Hd and t constant. From this we obtain: 

= (m| + mf) + |/ip) Hi + + l/rp) H^ - 2BHuHd , (D.38) 

M 3 = —2t‘^{AtHu — XtfiHd ), (D.39) 

j = 2 (2«" + i/J + i/J) , (D.40) 

A = + 2//3 + (//2 _ i/2 _,2y ,il 41) 


Dne to the large Yukawa coupling of the top quark, we also have to check the lifetime of 
the electroweak vacuum for sizable values of k, as opposed to the previous cases. Hence we 
can no longer approximate S{k) by *F(0), and we must minimize the following (normalized) 
bounce action: 


5 


A I /(^)-/(o) \ 

V ^(0) ) 


z^Ml 

~Mf^ 


(D.42) 


with [74] 


vrVe 16.5 28 

(1 — Ak)'^ (1 — 4k)‘^ 1 — 4/a: ’ 


(D.43) 


and 5(0) ~ 45.4. The minimization goes again over two variables, for example the ratios 
Hu/t and Hd/t, and is done numerically. It is important to keep in mind that S should 
be minimized only over the range 0 < /a: < 1/4, for which the potential is metastable. For 
/a: < 0 the potential is unbounded from below and for M| < 0 it is unstable, while for 
/a: > 1/4 the electroweak vacuum is the global minimum of the model. 


D.4 How to improve the estimate 

We can improve the estimates described above in several ways: 

• on top of the 6 cases described in Table 5, one could consider the case where both 0 
and 0 are leptons. One would then need to perform a minimization with respect to 
three helds, since there is no constraint due to SU(3) D-terms; 

• one could consider more general radius-dependent directions in held space, which 
would require the use of numerical algorithms as in Ref. [69]; 

• we have calculated transitions between vacua of the tree-level scalar potential at the 
hxed scale rrisusy To improve our estimates higher order corrections and/or RG- 
improved potentials should be considered. Such corrections can be relevant in some 
cases [61]. 

All these generalizations are however quite involved and are beyond the scope of this paper. 
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